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Chapter    I 
INTRODUCTION 

The  program  of  this  report  is  based  on  class  work  per- 
formed during  the  Summer  Quarter  1982  at  the  Naval  Postgrad- 
uate School  and  taught  by  Dr.  Stephen  W.  Tsai  of  the  Air 
Force    Wright    Aeronautical    Laboratories.  The   original    pro- 

gram was  based  on  reference  1  and  was  written  in  order  to 
compute  and  observe  a  variety  of  in-plane  and  flexural  par- 
ameters calculated  for  symmetric  laminates.  Upon  completion 
of  that  program,  however,  it  was  enhanced  by  the  addition  of 
subroutines  and  displays  to  compute  the  max  stress  and 
strain  strength  ratios  in  addition  to  quadratic  strength 
ratios. 

Chapter  II  explains  program  logic  based  en  the  flow 
charts  of  Appendix  A.  The  reader  will  find  that  no  special 
programming  techniques  were  required  and  all  calculations 
are    straight    forward. 

There  are  •'•wo  basic  advantages  to  having  a  program  such 
as  this  on  a  microcompter.  First,  although  a  microcomputer 
such  as  the  HP-S7  is  obviously  net  as  portable  as  a  hand- 
held calculator,  the  user  will  note  distinct  advantages  to 
being   able   to   view    labeled    output   and    matrices    several    lines 
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at  a  time  as  opposed  to  single  line,  unlabeled  output.  Sec- 
ondly, the  program  code  is  very  easy  to  read  and  modify  if 
necessary. 

This  program  is  nor  intended  to  be  all  encompassing,  but 
it  should  be  of  assistance  to  users  designing  symmetric  la- 
minates. Users  are  encouraged  to  modify  lhs  program  to  suit 
their  own  needs  and  to  publish  addendum  and  erratta  *o  this 
report.  Corrections  and  modifications  should  be  addressed 
to  Dr.  Tsai  at  AFWAL/MLBM,  Wright-Patterson  AF3,  Ohio. 


Chapter    II 
PROGRAM    LOGIC 

Ever,  though  this  program  has  been  written  fcr  the  HP-87, 
modifications  to  fit  this  program  en  other  microcomputers 
should  not  be  difficult.  Tha  only  really  unique  aspect  of 
the  program  occurs  during  the  initial  parts  of  the  program. 
Refer   to    the    flowcharts    of    Appendix   A. 

A  basic  familiarity  with  the  HP-87  is  assumed.  Load  the 
program  from  disk  storage,  ^hen  initiate  program  operation 
with  the  HUN  key.  The  program  will  prompt  the  user  to  sig- 
nify if  printer  output  is  desired  or  not.  Reply  with  Y  or 
N.  The  program  code  is  set  for  a  printer  address  of  791. 
Ensure  the  printer  has  this  aidress,  or  list  the  program 
code   and    change    the   address   in   line   60. 

The  HP-87  allows  a  programmer  to  set  up  14  predefined 
keys  for  selective  branching.  In  this  program,  selections 
are  established  for  each  of  five  commonly  used  laminate  ma- 
terials: T300/5208,  3(4)/5505,  AS/3501,  SP1002,  and 
KV49/EP0XY.  Selection  of  any  one  of  these  five  automatical- 
ly inserts  engineering  constants  into  appropriate  variable 
names.  The  units  of  these  variables  ar=  in  the  SI  system. 
Users   desiring   to    work    in    English    units    need    only    select    the 
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"SI  ->  SNG"  key  prior  to  selecting  a  material.  A  "flag  va- 
riable" is  set  to  unity  for  later  logic  testing  to  determine 
whether  or  not  to  apply  a  conversion  factor  to  each  of  the 
input  variables. 

Users  desiring  to  define  their  own  material  may  lo  sc  by 
selecting  "User  Def'd".  The  program  will  prompt  for  a  ma- 
terial name  and  each  of  the  engineering  constants  before 
continuing . 

The  "Reference"  key  selects  a  display  with  an  abbreviated 
explanation  of  the  program,  a  list  of  the  references,  and 
author  information.  Press  CONTINUE  to  return  to  the  materi- 
al selection  display. 

Unit  ply  thickness,  ho,  is  then  requested.  Ensure  that 
the  dimension  of  this  parameter  is  the  same  as  that  of  the 
orevious  variables. 

Following  this,  and  in  several  subsequent  locations,  the 
user  can  choose  to  display  computed  variables  or  to  bypass 
the  display.  If  a  printer  is  used  with  the  program,  prompts 
for  variables  and  displays  will  not  appear  on  the  CRT,  but 
will  appear  on  the  printer  instead.  Refer  back  and  forth 
between  the  two  to  determine  what  the  next  input  should  be. 
With  a  printer  connected,  only  a  question  mark  will  appear 
on  the  CRT.   First  time  users  are  encouraged  to  use  only  the 


CRT  display  (no  printer)  in  order  to  g=>t  used  to  answering 
prompts  and  observing  output.  Run  time  for  a  complete  ser- 
ies of  calculations  requires  only  a  few  minutes.  Errors 
should  be  handled  by  shopping  program  operation  with  the 
PAUSE  key,  then  reinitialize  with  the  RUN  key.  Effort  has 
been  expended  to  make  this  program  as  user  "friendly"  as 
possible.   Most  prompts  are  self  explanatory. 

Following  optional  displays  of  modulus  and  compliance  ma- 
trices as  well  as  strength  parameters,  tha  program  will  sig- 
nify that  computation  of  in-plane  variables  are  now  pro- 
gressing. A  prompt  for  CORE/NO  CORE  laminates  is  then 
displayed.  CORE  laminates  are  those  comprised  of  a  honey- 
comb sandwich  section  placed  on  the  centerline  that  contri- 
butes little  *o  the  stiffness  of  the  laminate  but  is  used' to 
alter  the  weight  characteristics.  The  thickness  of  the  cere 
is  usually  input  as  an  integral  multiple  of  the  ply  thick- 
ness.1 Proper  input  of  ply  information  in  sequential  order 
is  very  important  for  flexurai  properties  later  in  the  pro- 
gram. CORE  and  NO  CORE  materials  both  require  ply  informa- 
tion to  be  entered  from  the  center  ply  outward  in  consecu- 
tive  order.     No   attempt   has   been   made   to   group   ply 


1    A    non- 
tere 


on-integral  multiple  for  core  thickness,  r.c,  ma?  be  en- 
ed,  bu_  '■hs  initial  ply  index  number,  ti ,  will  still 
have  to  be  entered  as  an  integer.  Use  the  nearest  multi- 
ple rounded  off  to  the  nexf  highest  number.  Slight  errors 
may  occur.  but  the-/  should  be  significant  only  for  lami- 
nates   with    a    relatively    thick    unit    ply   thickness. 


orientations,  and  the  user  will  be  required  to  enter  each 
ply.  Programming  was  simplifies  by  this  expediency,  but  us- 
ers entering  complicated  laminates  or  laminates  with  many 
plies  may  find  it  somewhat  tedious.  Feedback  from  the  prin- 
ter is  immediate,  however,  so  any  mistakes  may  be  discovered 
quickly. 

Following  input  of  CORE/NO  CORF  information,  -"-he  in-plane 
modulus  of  laminates,  Rij,  is  computed  followed  by  its  in- 
verse, aij.  Normalized  versions  of  these  matrices,  Aij*  and 
aij*,  are  computed  as  well  as  the  effective  in-plane  engi- 
neering constants.  A  display  of  these  variables  is  option- 
al. 

The  program  then  requests  input  of  loads  (stress  resul- 
tants) N1,  N2,  and  N6,  followed  by  a  request  for  the  orien- 
tation of  the  ply  to  be  examined.  Variables  computed  from 
these  inputs  include  the  in-plane  strain  variables,  or.-axis 
strain  and  srress  variables,  max  stress-strain  strength  rat- 
ios (E/R ' ) ,  and  quadratic  strength  ratios  (R/R1).  k  non-op- 
tional display  lists  results  of  these  calculations.  Follow- 
ing the  display,  the  user  can  enter  a  new  load  and 
orientation  and  repeat  the  calculations  and  display.  If  a 
new  load/orientation  is  not  desired,  the  program  proceeds  to 
flexurai    computations. 


Without  further  prompts,  the  program  computes  the  flexur- 
al  modulus  matrix,  Dij,  followed  by  its  inverse,  dij.  And 
in  a  manner  similar  to  in-plane  computations,  normalized 
parameters,  Dij*  and  dij*,  are  computed  as  well  as  effective 
flexural  moduli,  Eif.  A  display  of  these  parameters  is  op- 
tional. 

The  program  now  prompts  for  the  moments,  M1#  M2,  and  [16. 
An  index  number,  t,  representing  the  sequential  ply  number 
for  the  ply  to  be  examined  (counting  from  the  ply  closes-  to 
"he  center  line)  is  requested  along  with  the  orientation  of 
that  ply.  With  these  inputs,  the  program  calculates  curva- 
ture, ki ,  on-axis  stress-strain,  max  stress-strain  strength 
rabies  (R  5  F')r  and  quadratic  strength  ratios  (?.  S  R')«  A 
non-optional  display  presents  all  parameters  in  a  format  si- 
milar  to    that    for    in-plane    parameters. 

Following  the  display,  the  user  can  enter  ?  new  set  of 
moments,  index  t,  and  THETAt  for  a  repeat  of  the  calcula- 
tions. Selecting  the  other  option  terminates  program  opera- 
tion . 


Chapter  III 
SAMPLE  CALCULATIONS 

The  following  pages  are  reproduced  from  sample  calcula- 
tions generated  by  the  LAMINATE  program.  The  first  program 
run  displays  all  optional  displays.  Subsequent  program  runs 
delete  many  of  these  displays  since  T300/5208  was  ased  for 
all  calculations  and  many  of  fhese  displays  would  not  vary. 
The  ply  designations  for  each  composite  are  noted  next  to 
the  title  fcr  each  run.  The  printer  used  for  these  runs  was 
configured  fcr  a  laterally  compressed  moie  to  allow  for  re- 
port format. 

Many  of  the  prompts  do  not  appear  on  these  print  outs. 
The  program  is  written  to  display  many  of  the  prompts  only 
on  the  CRT.  The  reader  is  referred  to  ths  flowcharts  to  re- 
solve questions  about  program  operation  a  specific  points. 


COMPUTATION  OF  IN-PLANE  AND  FLEXURAL  PROPERTIES  OF  SYMMETRIC 

LAMINATES  ON  THE  HP-87  rtlCROCOMPUTER 


Selected  conposite  natenal  =>  T3Q0/S208 
Input  laninate  thickness,  ho  (SI  *r  EN€), 

ho  =  ,000125 


[0/9  0} 


inmtmmtmttmuttttmKutittixtitutttttmxtxtttumttmnmttm 

Engineering  constants  and  compressive,  tensile  -  and  snear 
constants  for  T300/S208 

Ex  =  i8i.Q80E+»0?' 

Ev  =  i0.300E+009 

n  -  ,28 

Es  =  7.178E+009 

X  =  1508.M0E+0*: 

V  =  1SOOJ00E+006 
f  =      40.30IE+00« 

V  -    246.300E+006 
S  =      6S.900E+306 
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mmmtmmtmtmmtmitttmttmtmtmmummttwmtmnmt 

On-axis  Qij's,  Sij's,  UiQ's,  ftij's,  and  lij's  for  T380/5298 
(Note:  ij=xys), 


(  18i.8iE+0Q9  2.8969E+Q09    8 

on-axis  Qi j  =  (  18.346E+Q09    0 

(  7.1700E+889 

(  S.5249E-012  -i,547E-Qi2    8 

on-axis  3ij  =  (  97, 087F-012    8 

(  139.47E-812 

UiQ  =  76.3682E+009 
U2Q  =  35.732SE+009 
U3Q  =  19.71Q4E+0B9 

U4Q  =  22.6G74E+089 
USQ  =  26.8804E+009 

(  22.726E+QQ6  ,362i2E+Q96  0 

on-axis  Ai i  =  (  1.2933E+006  0 

(  .8962SE+3Q& 

<  44.199E-089  -12.37E-009    8 

on-axis  ai j  =  (  776.70E-039    0 

(  ii29,8E-899 

tmmttmmmttttmmtmmtmmmmtmumtmnnmtmttmm 


ttmttntmttmttutttttttttttmuittmttittmmttutitttuttttmtmttt 

Strength  parameters  in  stress  space,  Fij,  and  strain  soace,  Bij, 

for  T3G8/52Q8  (Note:  ij=xys), 


(   ,44444E-0i8    -3.360E-518  0 

Fij  =  (  ili.63E-0i8  0 

(  2ifc,26E-8i8 

Fx  =  B8J88E-801 

Fv  =  29.93SE-889 


(  12.B84E+0B3    -3.869E+883  8 

Bij  =  (  it,  saiE+003  0 

(  ii.ii8E+803 

Gx  =  6.0M7E+M1 
Gy  =  2.166QE+Q82 

imximmxmmtxtmxtmnmtmmmmtmnmmnmunnmmitt 
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mxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxtxtxxxxxxxxxxxxxxxxxxxxxxxxxx 
mxxxxxxxxxxxxxixxxxxxtxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxtxxxxx 

The  progran  will  now  conpute  in-plane  and  flexurai  prooerties  of 
laminates  specified  by  the  user.  Either  CORE  or  HO  CORE  lanmates 
can  be  specified  and  any  number  of  plies  defined.  The  user  will 
note  that  this  progran  takes  a  'brute  force'  approach  bv  reauest- 
mg  an  entry  for  each  ply.  This  method  allows  sinpler  proaran 
code,  ensures  fewer  errors,  and  enhances  accuracy. 

ttttuxxxtxxtxmxxxxxxxxxmxxxtxtxttxmtmxxxxxxxxxxxxxxmxxixxxxxxtxxtxtti 

txxxxxxxxxxxxxxxxxxxxxxxxxxxxxxmxxxxxxxxxxxxxxxxxxxxxxxxxxxxxixxxnxxxxxtxxxxx 

In-plane  coHoutations  follow.,. 

Hill  the  lawinate  contain  CORE  or  NO  CORE? 
C/NC  =)  END  LINE 


HO  CORE  laminate:  Enter  V,  the  total  nunoer 
of  plies,  surface  to  surface, 
n  =  ? 


NO  CORE  lonmate  thickness,  h  =  .0005 


Now  enter  the  piv  orientation  of  each  ply: 

THETA(i),  THETA<2),  .,,  ,  THETA(n/2) 

XXX  -->  FROM  CENTER  TO  OUTER  SURFACE  <-  XxX    (IMPORTANT!) 


THETA(  i  )  =  ? 

?0 
THETA(  2  )  =  ? 

0 
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mtxtttmtmmmtm&mnttmmtmttttuittmtmutmttmmtuut 

In-plane  paraneters:    ftij's,  dij's,  engineering  constants,  and 
normalized  Ai it's  and  aijt's  for  T300/52Q8    (Note:     ij=i26) 


(  48.Q39E+086  1.4485E+I06    8.4343E-906 

in-plane  Aij  =  (  48.039E+006    2.0144E-894 

(  3.535SE+006 

(  23.835E-00?  -.6282E-009    -,137iE-0i9 

in-plane  aij  =  (  29.83SE-009   -.1169E-017 

(  2.7394E-007 

Eio  =  9.599iE+0i0 
E2o  =  9,599iE+010 
v2io  =  ,03015 
E&o  =  7.1708E+009 

(  9.6079E+010  2.3969E+009    i . &869E-0Q2 

in-plane  Aij*  =  (  9.6079E+010    4J238£-00i 

(  7.1700E+009 

(  1.04i8E-91i  -.3i4iE-0i2    -.6859E-923 

in-plane  aij*  =  <  iJ4i8E-9ii    -.5846E-921 

(  i .  3947E-910 

tmttmtmttmtmttmtnttitutmmttmmmmmtmtmttmmmt 


Input  loads  rtl,M2,N6 


lii  =  i   N2  =  0   M6  =  9 


Mow ,  enter  THETAt,  the  orientation  of  the  ply 
*o  De  examined  (in  degrees), 
THETAt  =  ? 


THETAt  =  0 
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numtttutmmmmmtttmmmtutttnmtmttmmttmmmmm 

in-plane  strain:  eio  =  20.835E-9W 
e2o  =  -.6282E-009 
e6o  =  -.1371E-019 

For  THETAt  =0     ..,,     (T380/S2Q8) 

on-axis  strain:  eio  =  2.083SE-998 
evo  =  -6.282E-4i0 
eso  =  -,i371E-4i? 

on-axis  stress:  sxo  =  37.863E+902 
svo  =  S3.858E+9Q0 
sso  =-.9836E-010 


strength  ratios  8  V 

mx  strain  3.98E+08S  3.98E+99S 
nai  stress  3.96E+095  3.96E+085 
quadratic         3.41E+90S     5.54E+805 

average  stresses:  sto  =  6.3i88E+9!8 
sti'  =  1.1077E+009 
(=f (quadratic  R  i  R»)) 

ntxntitntxuiunuixtnitiiuxtnxmixtmuixixtutiinxtnumxuuiwt 


Hew  load;    Mi  =  i       f*2  =  11        M6 


How,  enter  THETAt,  tne  orientation  of  the  ply 
to  be  examined  (in  degrees), 
Tl 


THETAt  =  99 
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mxxxxxxxxxtxxxxxxxxxxxxxxxxtuxtxxxxxxxxxxxxxxxxtxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

in-plane  strain:  eio  =  20.835E-00? 
e2o  =  -.6282E-0Q9 
e6o  =  -. 1371E-019 

For  THETAt  =90     ,...     (T3M/52Q8) 

on-axis  strain:  exo  =  -.&232E-IJ99 
evo  =  2G.33S-G89 
eso  =  -,i?64E-»iB 

on-axis  stress:  sxo  =  -5.385E+0t)i 
syo  =  21.374E+091 
sso  =  -.i408E-Q08 


strength  ratios  R  R' 

nax  strain  L86E+0I5  i.iSE+Qufc 
wax  stress  i.87E+005  i.iSE+006 
quadratic         i.87E+QI5     1.13E+036 

averaae  stresses:  sto  =  3.7340E+908 
5t«>  =  2.2fc88E+089 
(=f (quadratic  R  &  R>)) 

ttmttmttttmttmtmttmmnmttttttmmtmtmmttmtnttmtmt 


mtmmmttmmnttmummntmtttmtmmttmtmttmtmtmm 

ixnxxxmxxtxnxxxtxxxmxttxxxxtxtxxxxxxxxxxxxxxxxxxxxxxxxxxxxxmxxxtxxnxm 

The  progra*  will  now  cowpute  flexoral  properties  for  the  ianmate 
defined  during  the  in-plane  definitions,    Properties  calculated 
include  the  flexurai  stiffness  Matrix  (Dij)»  the  flexural  compli- 
ance matrix  (dij)>  effective  flexoral  Ianmate  Moduli  (Eif),  and 
the  normalized  versions  of  stiffness  matrices  (Dijt  and  dij*), 

mixxxxxxxxxtxxxxxxxtxmxxmxnxxxxxxxxxxxxxxtxtxxxmxixxxxxtmxxxntumx 
txxxxxxmxxxxxxxxxxxxxxxxxuxxxxxxxxxxxmxtxxxxmxmxxxxxmxtmxxxxxxnm 
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tttmmmmmmnmtmttmmmmmmtmtmtnmtmttmttmxt 

Flexoral  properties  Dij,  dij,  effective  flexural  engineering 
constants,  and  the  normalized  Dijl  and  dij*.     (Mote:  ij=i26) 

(   1.670&E+080     .30176E-991    4.3929E-Q14 

flexural  Dij  =  (  3,3i03E-00i    i.4492E-Qi2 

(  7.4688E-002 

(  5,9957E-0Gi    -,5465E-O0i    4.iS12E-9i3 

flexural  dij  =  (  3.02S3E+000    -.4247E-010 

(  1.3389E+M1 

Eif  =  l.faOilE+Jli 
E2f  =    3.1727E+010 
v2if  =      .09116 
E6f  =  7.i705E+009 

(  1.6038E+011    2.3969E+089    4.2171E-003 
flexural  Dij*  =  (  3.i779E+0iQ    1 . Q072E-001 

(  7,l700E+t)09 

(  6.2456E-012    -.5c93E-0i2    4.3241E-024 

flexural  dij*  =  (  3.i5i9E-8ii    -.4424E-S21 

(  1.3947E-910 

%%mmmtxiitmmnmminmt%mmim%immmmm%tmttxtxmtt 


Input  nonents  ttl,M2>H6 


Hi  =  i       H2  =  0       *6  =  0 


Wow,  enter  the  index  nunber,  t,  and  *he  orientation,  THETAt 
)f  the  piy  to  be  examined. 
t,  THETAt  =  ? 


t  =  2   THETAt  =  I 
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mmttutttmmn%mttmt%xmmttmmx%ittimtnttimtutt%mMtn 

curvature,  ki  =  5 .99S7E-001 
k2  =  -.5465E-001 
k6  =  4.iSi2£-0i3 

For  t  =  2  and  THETAt  =  I 

flex'l  strains,  elf  =  i,4989E-G34 
«2f  =  -. 1366E-004 
ibi  '  1.0378E-016 

flex'l  stresses,  sif  =  2.7213E+007 
s2f  =  2.9286E+905 
sif  =  7.44iJE-007 

strength  ratios  R  V 

mi  strain  S.53E+001  5.53E+9Q1 
(tax  stress  5,5iE+00i  5.5iE+O0i 
quadratic         4.9SE+801      7. iiE+OOi 

averaae  stresses:  sto  =  i , 1880E+0I9 
sto'  =  i.705&E+Q09 
(=f (quadratic  R  U'i) 

mttummtttniuuttmmmtmmttimmtxmtmmttttumttmmt 


Input  new  loads  Hi,M2,H6 


Uew  iionents:  rti  =  i   H2  =  0   H6  =  0 


Uow,  enter  the  index  nanoer,  t,  and  the  orientation,  THETAt, 
of  the  ply  to  be  exanined, 
T,  THETAt  =  1 


x   =  t   THETAt  =  90 
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mmitmmmmwnmmmnmmittmnmtmtmmmtmtxmmi 

curvature ,  ki  =  5.99S7E-001 
k2  =  -.5465E-001 
k6  =  4. I5i2E-0i3 

For  t  =  i  and  THETAt  -  90 

flex'l  strains,  eif  =  7.4947E-I0S 
e2f  =  - . 6832E-005 

ebf  =  5.i890E-8i7 

flex'l  stresses,  sif  =  7.S562E+0OS 
52f  =  -.1025E+007 
s6f  =  -.5132E-00S 

strength  ratios  R  R' 

nax  strain       5.18E+J01      3.19E+002 
nax  stress       5.29E+801      3.26E4-I0S* 

quadratic         5.22E+001      3J1E+GG2 

a«eraoe  stresses:  sto  =  1.2535E+939 
sto'  =  7,2i39E+0Q9 
(=f (Quadratic  R  4  R')> 

nxxxitxtxmxxxxxxxxxxtxutxnxxxtxxxxxxtxxtxxxmtxxxxxxxxxxxxtxmxtxxxxxtxxx 


End  of  onaran.    Thank  you, 
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COMPUTATION  OF  IN-PLANE  AND  FLEXURAL  PROPERTIES  OF  SYMMETRIC 
LAMINATES  ON  THE  HP-87  MICROCOMPUTER 

Selected  composite  Material  =)  T3Q0/5208  lOp/^Op /45p  /~45p\ 


Input  lamnate  thickness,  ho  (51  «r  ENG). 


ho  =  .000125 


nmmttmnttmttmmtnttuttttmtmtuummtmtmtnMmmttt 

tttmmmttmmttmmmmmnmmitmummtmmt-intmmin 

The  progran  will  now  compute  in-plane  and  flexural  prooerties  of 
laninates  specified  by  the  user,    Either  CORE  or  HO  CORE  lanmates 
can  be  specified  and  any  ntmber  of  plies  defined.     The  user  will 
note  that  this  program  takes  i  'brute  force'  approach  by  request- 
ing an  entry  for  each  ply,    This  fiethod  allows  simpler  progran 
code,  ensures  fewer  errors,  and  enhances  accuracy. 

mtmmmuntmmttmttmmmtmtmmuttmttnmmttttnmtw 
mttmmmmtmmmttnuuMmmmmmtttmtmnmitmmtm 

In-plane  computations  follow,,, 

Will  the  lacunate  contain  CORE  or  tfO  CORE? 
C/HC  =)  END  LINE 


HO  CORE  lamnate:    Enter  V,  the  total  number 
»f  plies,  surface  to  surface, 
n  =  ? 
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=  16 


NO  CORE  laninate  thickness,  n  =  ,002 


Now  enter  the  ply  orientation  of  each  plv; 

THETA(i),  THETA(2),   ,..  ,  THETA(n/2) 

*tt  --}  FROM  CENTER  TO  OUTER  SURFACE  <-  »t  (IMPORTANT!) 


THETA(  1  )  =  ? 

-45 
THETA<  2  )  =  ? 

-45 
THETA(  3  )  =  ? 

45 
THETAC  4  )  =  ? 

45 
THETA(  5)*! 

90 
THETAt  6  )  =  ? 

90 
THETA(  ?)»1 

THETAi  8  )  =  ? 
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mmttmmttumtmtmmmmttmtnmmtmttmmmttttmmttt 

In-plane  parameters:    Ai.j's,  aij's,  engineering  constants,  and 
nornalned  Aij*'s  and  aijt's  for  T300/5208    (Note:     ij=l26> 


(  1S.274E+907  4.S215E+0O7    1.68691-005 

m-olane  Aij  =  (  iS.274£+Q07    4J288E-0M 

(  S,376lE+907 

(  7i,76iE-0i9  -.2124E-908    i.3668E-02fl 

in-plane  aij  =  (  71.761E-010    -.5311E-01V 

(  i.3601E-008 

Eio  =  6.9676E+010 
E2o  =  6.9676E+M 
v2io  =  ,29653 
E6o  =  2.6880E+010 

(  7.6368E+01I  2.2607E+910    8.4343E-003 

in-plane  Aij*  =  (  7 . 6368E+0 i 0    2 , N44E-Q01 

(  2,6380E+8i9 

(  i.4352E-0ii  -,4248E-0ii    2.7334E-023 

in-plane  aij*  =  (  i,4352E-9ii    -,i062E-l2i 

(  3.7202E-811 

mtmmmttttmmmtwtmtttmmmmtmxtMttmtmtmmtttm 


Input  loads  Ni,N2,N6 


Mi  =  1       M2  =  I       N6  = 


Mow,  enter  THETAt,  the  orientation  of  the  ply 
to  be  exatuned  (in  degrees). 
THETAt  =  ? 


THETAt  =  45 
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mtmttmtmxtxttmmmmtmtttttntmtmtmitxtinmmtitmiitn 

in-plane  strain:  eio  =  7i .7&1E-918 
e2o  =  -,2i24E-308 
e6o  =  i  .3668E-920 

For  THETAt  =45    ,,..     (T300/S208) 

on-axis  strain:  exo  =  2.5259E-089 
en  =  2S.259E-019 
eso  =  -.9300E-003 

on-axis  stress:  sxo  =  46.t>55E+00t 
syo  =  33.4S9E*tOQ 
sso  =  -.6668E+002 


strength  ratios  R  R' 

«ax  strain  1.02E+006  1.02E+006 
fwx  stress  1.02E+006  i.82E+08b 
quadratic         &.94E+0Q5      i.35E+006 

average  stresses;  sto  =  3 . 4?98E^00a 
sti»  =  6.7538E+008 
(=f (quadratic  R  h  R»)) 

mnmtttmtmimnmtmtmttntnmmtumtmmnmmitmimti 


Hew  load;  Hi  =  i   N2  =  3   N6  = 


Xow,  enter  THETAt,  the  orientation  of  tne  ply 
to  be  examined  (in  degrees). 
THETAt  =  ? 


THETAt  =  '0 
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mtntmtmmmmttxtttmmmnmtmtmmtmmtmnmttutmtt 

in-plane  strain:  eio  =  7i .76iE-0it 

e2o  =  -.2124E-038 
efao  =  1.3668E-020 


For  THETAt  =  90 


[T3I8/S2081 


on-axis  strain',  exo  =  - . 2124E-0Q8 
eyo  =  7i,76iE-Qi9 
=so  =  -.1347E-018 

on-ixis  stress:  sxo  =  -3.65AE+082 
svo  =  :8J9iE+000 
sso  =  --7S09E-089 


strength  ratios 


R 


rt> 


mox  strain  5.41E+3QS  3.33E+006 
«ax  stress  5.37E+O0S  3.61E+006 
quadratic         S.52E+305      2.oOE+006 

averaae  stresses:  sto  =  2 . 7612E+088 
sto'  =  1.2976E+909 
(=f (quadratic  R  4  *')) 


^>w  Iqqo;    Mi  =1       HZ  -  I        N6 


Now,  enter  THETAt,  the  orientation  irf  the  ply 
to  bi  examined  (in  degrees). 
THETAt  =  ] 


THETAt  =  0 
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mtmxtmttxxtxxxxxxxxxttxtxxmxxtxtxxxxxxxxxmxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


in-plane  strain:  eio  =  7i.7ME-8ii 
e2o  =  -,2124E-!>08 
e6o  =  i.3668E-«2li 


For  THETAt  =  8 


(T3M/S208) 


on-axis  strain:  exo  =  7. 1761E-00? 
evo  =  -2.124E-009 
eso  =  1.3668E-I25 

on-axis  stress:  sxo  =  12 . 985E+C02 
svo  =  -i .  190E-»-OaO 
sso  =  9.7999E-011 


strength  ratios 


? 


? 


mqx  strain 
wax  stress 
quadratic 


i.i5E+906     i  .i3E+906 

1.16E+006      i.i6E*0«l. 
i .I6E+906     i .13E+006 


average  stresses:  sto  =  5.3l3iE+0S9 
sto'  =  5.6542E+QI8 
(=f (quadratic  R  &  R')) 

mmmtmmmmmttmttnmmnmmtmmmmtmmumtitxiu 


mxmttnttmttnnmimntmtmiutmmmxttmtxntmmmmtiut 
mxttxxxxxxxxxxxxxxxtxxxxxxxmxxxtxxxxxxxxxxnmxxxxxxxxxxxxxxxxxxxxxxnxxxxx 

The  qroqran  will  now  cowqote  tlexural  orooerties  for  the  lamnate 
defined  during  the  in-plane  definitions,    Prooerues  calculated 
include  the  flexurai  stiffness  natnx  (Dij),  the  flexorel  compli- 
ance natrix  (dij),  effective  flexurai  latunate  moduli  (Eif),  and 
the  normalized  versions  of  stiffness  Matrices  (Dij*  ana  dij*). 


xxxnxxxxxxxxxmxxxxxxxxxxxxtxxxtixmxxxxxxxxtxxxxxtxtxxxxxxxxtxxnxxxxxxxxxxx 

xxxxxxxxxxtmxxmxxxxxxxxxxxxxxxxxtxxxxxxxxxxtxxxxtxxMxxnxxiXimixxmmif. 
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mmmmtmtmmmttmmumttmtmmtmmummmmtttmtt 

Flexural  properties  Dij,  dij,  effective  flexural  engineering 
constants,  ana  the  nornalized  Dij*  and  dij*.     (Note:  1  j=126) 

(  7.6342E+091     .52iME+08i    2.6791E+000 

flexural  Dij  =  (  4.4693E+98i    2.679iE+Q00 

(  8.0651E+000 

(  i.3241E-9G2    -.1307E-9G2    -.39ME-902 

flexural  dij  =  <  2.2959E-002    -.7192E-002 

(  i,2779E-9Qi 

Elf  =  i.i329E+0ii 
E2f  =    6.5334E+010 
v2if  =      ,39877 
E6f  =  1.1747E+010 

(  i.i526E+9ii    7.8245E+M9    4.9137E+909 
flexural  Dij*  *  <  6.7B39E+0t0    4.8187E+009 

(  i.2898E+0i0 

(  3,827iE-9i2    -,87i8E-9i2    -,2&42E-9ii 
flexural  dij*  =  (  1.S396E-3U    -,4794E-9ti 

(  8,5i32E-9ii 

mmm**t*m*mtm*t*x*mm*m*mi***m**t*m*mm*t**m**t*m*s 


Input  nowents  Ki,H2,H6 


Hi  =  i       H2  =  9        H6  =  9 


U«w,  enter  the  index  nuwber,  t,  and  the  orientation,  THETAt, 
of  the  plv  fo  be  examined, 
t,  THETAt  =  ? 


t  =  3   THETAt  =  9 
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ttummntunuummtmmmtmmmitixtttmmuttmtntttmmt 

curvature,  ki  =  1.3241E-002 
k2  =  -. 1307E-002 
U  =  -. 3964E-062 

For  t  =  8  and  THETAt  =  I 

flex'l  strains,  eif  =  i .3241E-00S 
e2f  =  - .  1307E-Q0S 

ebf  =  -.3964E-905 

flex'l  stresses,  sif  =  2 . 403SE+006 
s2f  =  2.4827E+304 
sif  =  - . 2B42E+005 

strength  ratios  R  R' 

«ax  strain       6.26E+082      S . 2SE+092 
nax  stress       6.24E+002      6.24E+90c: 

quadratic         5.46E+Q02     7.&2E+0G2 

average  stresses:  sto  =  3,i38iE+Q8S 
sto'  =  i.i43i£+88? 
(=f (quadratic  R  4  R')) 

mxtixxmxtxxxxxxtxtxxxxxittixttmtxxxxxxtxtxtxtxtxxxxxxxxxtxxxxxxxxxxxxmxx 


Input  new  loads  rti  ,{42,116 


Mew  Nentsi    Hi  =  i       H2  =  I       ,1S  =  3 


^ow,  enter  the  index  nunber ,  t,  and  the  orientation,  THETAt 
of  the  oli  to  oe  exanined. 
t,  THETAt  =  ? 


t  =  6       THETAt  =  91 
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mtmmttnmtmm%tmmtttmm%txttmntx*ttnimtm%tttmmmt 

curvature,  ki  =  i.324iE-002 
k2  =  -. 1307E-002 
k6  =  -.3964E-002 

For  t  =  b  and  THETAt  =  90 

flex 4  strains,  eif  =  9.939SE-006 
e2f  =-.?808E-006 
eif  =  -.2973E-995 

flex'l  stresses,  sif  =  ?,??01E+004 
s2f  =  -.149SE+006 

s6f  =  -,2i3iE+00S 


strength  ratios 


R 


wax  strain 
«ax  stress 
Quadratic 


3.91E+002  2.41E+003 
4.00E+002  2.46E+Q03 
3.39E+0I2     2.10E+003 


average  stresses:  sto  =  5.3413E+099 
sto>  =  3.1496E+099 
(=f(quadratic  R  4  R')) 

nmummttttmmmttnmmtwttttmttmtmtmmttmttttxttuttt 


Input  new  loads  Mi.fl2,rt6 


^'ew  .^o«ents:    (11  =  i       rt2  =  i)       ?16  =  3 


Nsw,  enter  the  index  number,  t,  and  the  orientation,  THETAt ; 
»?  +he  Dly  to  be  examined. 

t  THETAt  =  ? 


t  =  4   THETAt  =  45 


27 


mttntmtmnmtmumnmmmmumumnmtmumtmumtm 

curvature,  kl  =  1.3241E-0Q2 
it2  =  -.1307E-092 
kS  =  -.3964E-002 

For  t  =  4  and  THETAt  =  45 

flex'l  strains,  eif  =  6,6203E-flfl6 
e2f  =  -.6538E-«0h 
e6f  =  -.1982E-405 

flex'l  stresses,  sif  =  2.6246E+005 

s2f  =  i,S315E+30S 
s6f  =  -.29S3E+006 


strenath  ratios 


3 


«ax  strain 
mx  stress 
quadratic 


9.77E+002  1.30E+I03 
3.53E+0G2  i  .3IE+0Q3 
6.71E+002      i .77E+003 


average  stresses:  sto  =  i . 0072E+089 
sto'  =  2.9548E+809 
(=f< Quadratic  R  h  V)) 

mtmmitnmttmtmmtttmmntititmuutmtMmmmtmtmtu 


End  of  proqran,     Thank  you. 
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COMPUTATION  OF  IN-PLANE  AND  FLEXURAL  PROPERTIES  OF  SYMMETRIC 
LAMINATES  ON  THE  HP-37  MICROCOMPUTER 


02/902/452/-452/Ce 


Selected  composite  mterial  =)  T30O/S2O8 


Input  iawinate  thickness,  ho  (31  ir  E'NG). 


ho  =  .000I2S 


mtnmtmtmtmmttmtmmmmtttumnmmtmmmttmtmm 
mtntmmmumttttmummtnmtmmutmumttntmnmMtttt 

The  proqraw  will  now  co&pute  in-plane  and  flexural  prooerties  of 
laeunates  specified  by  the  user.     Either  CORE  or  NO  CORE  laminates 
can  be  specified  and  any  nunoer  of  plies  defined.     The  user  will 
note  that  this  program  takes  a  'brute  force'  approach  ov  request- 
ing an  entry  for  each  ply.    This  Method  allows  sinpler  proqran 
code,  ensures  fewer  errors,  and  ennances  accuracy, 

mtmtmmttunttttmnttmmtttmmttnmmtmttmmtmmtm* 
mtttxmmttmtmmmmttmttxmmummmmtmntttmmmm 


J  s 


In-plane  computations  follow... 

Mill  the  Lamnate  contain  CORE  or  NO  CORE? 
C/'NC  =)  END  LINE 


CORE  lamnate :     Enttr  V,  the  total  nunber  of 
lamnate  plies,  and  'nc'.  the  total  nuaber  of  equivalent  core 
plies  that  define  its  thickness,  surface  to  surface. 
n,nc  =  \ 


16 


nc  =  16 


CORE  hamate  thickness,  h  =  ,304 
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Ciuntinq  outward  fro«  the  center  of  the  latunate,  enter  an 
integer  number  'ti'  (t  initial)  that  represents  the  consecutive 
index  nunber  of  the  ply  adjacent  to  the  core.     (For  instance, 
if  the  core  is  a  total  of  4ho  thicx,  count  2  equualent  ply  layers 
fron  the  center  so  that  ti=3). 
ti  =  ? 


ti  --  9 

H«w  enter  the  ply  orientation,  in  decrees,  for  each  ply: 

THETA(ti),  THETA(ti+i)  ....  THETA(n/2>. 

*«  -)  CENTER  MOST  PLY  TO  OUTER  SURFACE  PLY  (-  *«     (IMPORTANT) 

THETA(  ?  )  -  \ 

-45 
THE7A(  it  )  =  ? 

-45 

THETA(  ii  )  =  ? 

45 
THETA(  12  )  =  7 

THETA(  13  )  =  ? 

90 
THETA(  14  )  =  ? 

?0 

THETA(  15  )  =  ? 

3 
THETA(  16  )  =  ? 

0 
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mmmmtmttittmtttmmtwmmnmmmmmttmmttmmmt 

In-plane  paraneters:    Aij's,  aij's,  engineering  constants,  and 
normalized  ftij*'s  and  aij«'s  for  T308/S238    (Note:     i j=i26) 


(  1S.274E+997  4.S21SE+497    1.6369E-O0S 

in-plane  Rij  =  I  15.274E+907    4J288E-904 

(  5.37ME+0Q7 

(  71.761E-010  -.2124E-998    1.3668E-920 

in-plane  aij  =  (  7i.76iE-010    -.S3iiE-li5> 

(  1.8601E-908 

Eio  =  3.4838E+010 
E2o  =  3.4838E+010 
v2io  =  ,29693 
E6i  =  i,3448E+810 

(  3.8i84£+0i0  l.i304E+010  4.2171E-003 

m-plane  ftijt  =  (  3.8i34E+OiO  1 . 8072E-901 

(  i.3440E+8i9 

(  2.3784E-0ii  -,84?7E-9ii    5.4672E-823 

m-plane  aijl  =  (  2.8794E-0il    -,2i24E-32i 

(  7.4484E-81I. 

tmmtttmtmntttumttmttmxmtmnmmttMmmmmmttmtt 


Input  loads  Ml >N2>N& 


Hi  =  i   N.2  =  )  %  =  9 


Now-  inter  TKETAt,  the  orientation  of  the  ply 
to  be  examined  (in  degrees). 
THETAt  =  ? 


THETAt  =  45 


txxxmxtxxxxxxtxxxxxxtxtxxxixtxmsxxxxxxxxxxxxxxxxxixxxxxxxxxxxxuxxxxxxxxxxtx 


in-plane  strain:  eio  =  7i.76iE-0iJ 
e2o  =  -.2124E-008 

i6o  =  i.36b8£-020 


For  THETAt  =  45 


(T338/S2Q8) 


on-axis  strain;  exo  =  2.S259E-00? 
evo  =  2S.2S9E-010 
eso  =  -.9308E-008 

on-axis  stress:  sxo  =  46.6SSE+051 
590  =  33.450E+000 
5S0  =  -.6668E+002 


strenath  ratios 


R 


i 


mx  strain  1.32E+006  1.02E+086 
«ax  stress  1.021+006  1.02E+086 
quadratic         6.94E+085      1.3SE+006 

averaae  stresses:  stu  =  i . 73SBE+0S8 
Sto'  =  3.37S4E+908 
(=f (quadratic  R  &  R')) 

ntmummimxtxtmttmtmmnmtmmimmmmxmntmttmm 


mxxnxtixxxxxxtmtixxxxxxmxttxxxxxxxxxxxxxxxxxmxxxxxxxxxxxxxxxxtxxixxxxxx 
xxxxxxxxxxxxxxixxxxtxxxxxxtxxxxxxtxxixxxxxmxxxxxxxxxxxxxxxxxxxxsxxmxxxxxxtxx 


The  orooran  will  now  compute  flexural  properties  for  the  lamnats 
defined  during  the  in-plane  definitions,    Prooerties  calculated 
include  the  flexural  stiffness  natrix  (Dij),  the  flexural  compli- 
ance natrix  (dij),  effective  flexural  larunate  Moduli  (Eif),  ind 
the  normalized  versions  of  stiffness  matrices  (Dij*  and  dij*). 

fxxxxxx*xtxxtmx*xx*xxxxxxxxxxx*xmtxxxxxxxxxmxxxtx*xmxixtxxxmxxxxxxxxxx 
tmxxxixxxxxxxxxxxxxxixxxxttxxxxxxxmxxxxxxxtsxxxxixixxxxxxxxxxxxmiixxmxxx 
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ttmmmmumtttmmtmtmmmmtmttummtttixmmmtmttt 

Flexural  properties  Dij,  dij,  effective  flexural  engineering 
constants,  and  the  nornained  Dij*  and  dijl.     (Note;  ij=i26) 

(  4.2346E+002    .7S935E+9Q2    1.3396E+081 

flexural  Dij  =  (  3.4844E+002    1.3396E+001 

(  9.S376E+001 

(  2.4645E-O03    -.S2&E-903    -.2797E-093 

flexural  dij  =  (  2.9930E-003    -.34S2E-003 

(  i,8Sit>E-G02 

Eif  =  7.6880Et0i9 

E2f  =  6.2S42E+010 

v2if  =  .21370 

E6f  =  i.7830E+0i3 

(  7.9398E+013    i .4238E+010    2.5ii7E+009 
flexural  Dij*  =  (  6.5333E+010    2.5117E+089 

(  i,7?77E+0i8 

(  i,3i44E-0ii    ~.2808E-Bii    -,i444E-0ii 

Flexural  dij*  =  (  i,5989E-0ii    -.i84iE-lii 

(  5,688bE-lii 

mtmtmtttmmittmnnmnmtttmmtmutmmmttummttxtttt 


Input  ^onents  Ni,N2,H6 


hi  --  1       H2  =  0       M6  =  0 


Miw,  enter  the  index  mnoer,  t,  and  the  orientation,  THETAt, 
»f  ^e  ply  to  be  examined, 


t,  THETAt  = 


f  =  16       THETAt  =  0 
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mtmtxmmtmmtnmttmmtmmttwnmtmxntmmtminntm 

curvature,  kl  =  2.464SE-903 
k2  =  -.5266E-0C3 
k6  =  -.2707E-003 

For  t  =  16  and  THETAt  =  I 

Plex'1  strains,  eif  =  4.9290E-006 
e2f  =-.i653E-085 

e6f  =  -.5415E-0Q4 

flex-1  stresses,  sif  =  8.9315E+00S 
s2f  =  3.3809E+003 
s6f  =  -.3882E+004 

strength  ratios  8  R> 

«ax  strain  1.&8E+903  i .b3E+003 
nax  stress  1.68E+093  i.o8E+Q83 
quadratic         i,o2E-*-003      1.83E+0I3 

average  stresses:  sto  =  4. 0642E+0I6 
sto'  =  6.3480E+008 
(=f (quadratic  8  4  8')) 

mmmimtmnimnmtiimmnmutttmunntttimtmtmmmtit 


Input  new  loads  !li,tt2,fffc 


d«w  iiwents:    Hi  =  i        H2  =  8        H6  =  9 


Miw,  *nter  the  index  nunber,  t,  md  the  orientation,  THETAt: 
of  the  ply  to  be  examined . 
t,  THETAt  =  ? 


t  =  14       THETAt  =  90 
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mtttmxnmtmmmtnmtmmmttmmtMtmtmmmtmxmmtt 

curvature,  ki  =  2.4645E-003 
k2  =  -.5266E-903 
k6  =  -.2707E-083 

For  t  =  14  and  THETAt  =  9t 

flex'l  strains,  elf  =  4.3129E-806 
e2f  =  -.9216E-006 
eif  =  -.4738E-Q86 

flex'l  stresses,  sif  =  4, 19S2E+004 
s2f  =  -,i55!)E+906 
s6f  =  -.3397E+004 

strength  ratios  R  R' 

wax  strain       9.30E+0Q2     S.54E+003 

«ax  stress       9.S3E+082      S.36E+G03 
quadratic         9.14E+082      4,o4E+Q83 

average  stresses:  sto  =  3.4234E+088 
sto'  =  i.7398E+f)09 
(=f (quadratic  R  4  R>)) 

mmtmtmttmttmmttmnmmmmtnttmttmmtmtwxttttmtt 


Incur  new  loads  Mi,rt2>M6 


:4«w  nawents:    Ml  =  1       K2  =  ■)       M6  =  0 


4jw,  enter  the  index  number,  t,  and  the  orientation,  THETAt, 
of  the  ply  to  be  exanmed. 
t,  THETAt  =  ? 


12       THETAt  =  45 
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mtmtmttmmmttMtmmmmmtmmmtttmmumttmmmu 

curvature,  ki  =  2. 4645E-003 
k2  -  -.526oE-003 
k6  =-,27l7E-40i; 

For  t  =  12  and  THETAt  =  45 

flex'l  strains,  elf  =  3.6968E-006 
«2f  =  -.7908E-006 
e6f  =  -.4061E-006 

flex'l  stresses,  sif  =  1.5361E+08S 
s2f  =  9.4270E+0G4 
s6f  =-.7019E^O!l 


strenath  ratios 


R 


\- 


«ax  strain 
nax  stress 
quadratic 


2.UE+0O3  2.UE+003 
L.93E+Q03  2,iiE+0Q3 
1.30E+003     2.97E+803 


average  stresses:  sto  =  4.3649E+Q18 
sto'  =  1.1153E+M9 
(=f(quaaratic  R  4  V)) 

mxtmtm%tmm*%%tm%i%mmimmmtmmmmmnmmttmmn 


End  of  praqran,     Thank  jou, 


Chapter    IV 
CONCLUSIONS 

The  original  intent  of  this  program  was  to  duplicate  and 
enhance  the  TI-59  program  of  reference  1  on  a  popularly 
available    microcomputer.  Ths    HP-87      is    a      unit    more      than 

likely  to  find  utility  in  many  engineering  organizations  and 
it  was  the  hope  of  the  author  to  contribute  to  the  design  of 
laminates  by  automating  many  tedious  calculations,  present- 
ing   them    in   an   easy   to    read  format. 

The  program  listing  in  Appendix  B  should  allow  guick  and 
easy  trouble  shooting  of  any  uncorrected  programming  errors. 
Modifications  and  additions  should  be  relatively  simple  to 
program. 

Overall,  "-his  program  allows  considerably  more  flexibili- 
ty than  the  TI-59  program  It  should  also  serve  well  as  an 
educational  tool  for  students  and  engineers  desiring  to  re- 
view   the    thecry    of   symmetric    laminates. 
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Appendix   a 
PSOGR&a    FLOWCHARTS 


InfttaMze 

variable*. 

constants,  ana 

arrays 


Set  key  labels 

for   ■aterlal 

s«lect?on 

display 


Convert  SI 
constarts  to 

ENG  urtts 
(Z6894. 9/543) 


<D 
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ieer  I  ng\__ 
tant*J>^^ 


-  enqtneer 
^"--.fons 


OtsDiay  Ex.  Ey. 
v«,  ES.  X.  X*\ 


NO 


:o*oute  on-a*ls 

modulus  and 

conpl tance: 

Of J.  Sf J 


eqns  (1.9).  (1.12)  00.13. 15 


Comoute 

O?Q*s=f(0xy*) 

t-1-5 


eqn  (3.15).  D.73 


Compute  on-axfs 

A  matMx  for 

single  Dly 

Al J-Qf ;<mo 


Compute 
atJ-AtJ**-l 


Compute 
strangtn 

oara«»eters  FT  J 
»nd  Gt  J  (f J=xy») 
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L 


Enter  n.nc  (n*4 
piles;  nc-4 
equlv  core 
pi  tes) 


J   h* (n*nc)*no 


Enter  • I ,  T ide* 

#  for  ply 

adjacent  core 


^- 


.^BH  THET* 


enter 
...  T^TA^n/2J 


1 


zs«nc*no/h 


9 
L_ 

Computation  of 

'  n-  d  l  ane 

properties 

fol lows. . . 


CORP 


HO   CORE 


Enter  n    (n-#  of 
piles) 


I SuBR  V.VVfj) 

H V1*2*ho*Vl 
VfO»«M**3)*?/3*V| 


volume  fraction  of  core 


Comoute 
tn-D'ana 
s^odu'us  of 
lamfnetes,  aIJ  j 

(MM  26) |   T4.3.  p. 124 


Comoute 

In- p 1 ane 

comp l 'ance 

!( Inverse)  »>  al  J 


L 


matrt*  Inversion 
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Compute 
nornai  l2e<J 
varfaoies: 
At 1»»a! f/h 


Compute  effec- 
tive fn-plane 
engineering 
constant*:  ETo, 

Yf.fq 


O 


Input  loads: 
Ml,  N2.  N6 


Compute 
tn-o'ane  strain 

variables: 
etc-ftatJ.Nt) 


stew  load  and  the  I  At. 
from  next  page. 


Enter  IHETAt. 

angle  of  ply  to 

5e  examined 


display  Atj. 

at  J.  Eto,  AfJ., 

■  if* 


TC4.2).  p. 121 


Compute  on-axfs 

strain 
var tables,  eto 


T2.6.  p. 53 
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t 


Compute  on-exts 

stress 
var'aDles.   sTo 


T1.6,    p. 14 


SWOT  ■iS'iR 


Compute  «n 

stress-strafn 

strength 

ratfos.  Ufl' 


QUAQ 


compute 
quadratic 
strength 

rat 'as.,  3,  ft  B! 


Display 

tn-piane 

strain,  on-sxfs 

stress-stratn, 

Ufl' 


Co  to  previous  page 


Computation  of 

f lexural 

properties 

follows  ... 


t 


^S-(l-ZS**3)*r.». 
f       3/1? 


See  T5.3,  0.180 


Compute 

f lexural 

moou  1  i-s . 

1 J»f (Uf0.ns.vf0 


T5.3.  p.190 
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coi«out*  aft 
( tnverse  ot J) 


Compute 
tj«-df  >h«»3/i; 


Compute 

elective 

f lexurai 

laatnatt  moCui  t 

ETf , .  yT  tf 


Eqn  (5.21).  p. 177 


Y©« 


01  splay  DM, 

Etf.  orj*. 

dtj* 


dij. 


Input  tiomenta: 

Ml,     M2.     M6 


Enter    fnoex  t 
and  TMETAt   for 
"H       <»»fred   ply 


New  »oment* 
(see  next  page) 


Compute 

curvatures 
kf-f(dfJ.Mj) 


J   T5.2.   0.175 
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x 

Co*Dute  on-axts 

stress-strain, 

elo  and  sto 


f?.6.    0.53  T1.6.    D.14 


.Syftr  BUAfl 


iycr. 

compute 

quadratic 

strength 

ratios,  ha* 


Display 

curvatures, 

on-ax?s 

stress-strain, 

DM',    for 

t/THETAt 


^>  "  >Q 


Go  to  orevlous  page 


End  of 
program 


Appendix   B 
PROGRAM     CODE 


13  PAGESIZE  16 

23  OPTION  BASE  1 

30  CLEAR  8  DISP  *  DISP  "Do  vob  wish  to  hove  a  copy  of  the  output  on  ' 

40  DISP  "a  printer?  Y/N' 

SO  INPUT  A* 

60  IF  UPC$  (A$)="Y"  THEN  PRINTER  13  701  ELSE  PRINTER  IS  i 

71  CLEAR 

30  PRINT  ?  PRINT  "COMPUTATION  OF  IN-PLANE  AND  FLEXURAL  PROPERTIES  OF  SYMMETRIC" 

90  PRINT  '  LAMINATES  ON  THE  HP-87  MICROCOMPUTER"  §  PRINT  ?  PRINT 

100  DISP  8  DISP  i  DISP  &  DISP  "Initialization  »f  prograw  enables.    Standby.3 

110  !     Initialize  variables,  strings,  and  arrays. 

120  INTEGER  I,J,t,ti,Mc,FlAG 

130  REAL  EX)H)NUX,ES,X)XP,Y)YP,S,eio)e2fl,e68,eif,e2f,e6f)5if,s2f,56f 

140  REAL  r1,QXX>QYY>«XY,QSS?SXX>SYY,SXY,^S,Ei»,E2o>v2io?E6o,Eif,E2f,92ifJE6f 

ISO  REAL  AXX,AYY,AXT,ASS,Mi,at»,axy,ass 

168  REAL  FXX,m,FXY,FSS)FX,FYJFXYS,GXX,GYY>GXY,gSS,GX,GY,UlQ,U29,l)3Q,U4Q)US< 

170  REAL  Hi,N2,N6)Mi,M2,H&,ki,k2,k6,h,ho)z5,hs,p,q,a,lj,c,nti,HH,SE,Ie,Re,x 

130  REAL  Vi,V2,V3,V4,VlD,V2M3D, W,THETA,THETAd >DET 

190  HEAL  £X8,eyo,e5*,sxo,5i|8,S50,5t8,5topp 

200  REAL  REPSX,REPSXP,REPSY,REPSYP,REPSS,RE?SSP 

210  REAL  RSIGXjRSIGXPjRSIGY^SIGYPjRSIGS^SIGSP 

228  DIM  A(6,6),a(6,6),AS(6,6),as(6,6),Qt6,6),THE7A(S0) 

230  DIM  D(6,6),d<6,6) ,58(6,6) ,ds(6,6) 

240  Vl,V2,V3,V4,VlD,V2D,V3D,V4D,FLAG=li 

2S0  SE=6894 .97543 

260  FXYSs-.S 

273  FOR  1=1  TO  6 

230       FOR  J=i  TO  * 

290  9(I,J)=0  3  A(I,J)=0  9  AS(I,J)=0  %  a(I,J)=0  §  as(I,J)=0 

300  D(I,J)=8  9  DS(I,J)=8  9  d(I,J)=0  *  ds(I,J)=8 

310       NEXT  J 

329  NEXT  I 

330  GOTO  660 

340  CLEAR  *  FLAGM 

3S3  PRINT  8  PRINT  "SI  -)  English  Units  Conversion;    A  progran  ?flaa;  has  been* 

360  PRINT  'that  will  key  appropriate  logic  to  convert  each  of  the  values  inpit" 

370  PRINT  'aotortatically  in  SI  units  by  progran  READ  statements.    If  this3 

330  PRINT  'option  is  selected,  be  sore  to  input  ;'bo\  the  laninate  thickness, '' 

390  PRINT  "in  English  units  as  well." 

400  PRINT  ?  PRINT  'Now,  press  CONTINUE  to  return  to  the  CRT  display  «f  the" 

410  PRINT  Material  selection."  9  PRINT  8  PRINT 

420  PAUSE 


-    UH     - 
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430  GOTO  660 

440  PAGESIZE  24  9   CLEAR  9   PRINT  9   PRINT 

450  PRINT  '   The  LAMINATE  program  is  designed  to  quickly  calculate  in-plane" 

460  PRINT  'and  flexiral  strength  and  stiffness  properties  of  selected  " 

470  PRINT  'laminates.  The  program  is  interactive  and  prompts  the  jser  for" 

480  PRINT  "ail  required  input  parameters.  Other  prompts  will  request  the  user* 

490  PRINT  'to  select  or  skip  a  display  of  various  modulus  and  compliance  ' 

500  PRINT  'matrices  as  well  as  effective  engineering  constants.  Host  other* 

510  PRINT  'prompts  are  self  explanatory' 

520  PRINT 

530  PRINT  "References:     (a)  Tsai,  S.U.,  Hahn,  H.T.,   'Introduction  to  Composite' 

540  PRINT  "  Materials',  TECHNOHIC  Publishing  Co...  ci?8fl." 

550  PRINT  ■  (b)  Donaldson,  S.L.,  'Revised  Instructions  for  TI-59* 

560  PRINT  '  Combined  Card/Module  Calculations  for  In-Plane  and' 

570  PRINT  '  Flexural  Properties  of  Symmetric  Laminates',  Air8 

580  PRINT  *  Force  Systems  Command,  June  1982.' 

590  PRINT  §  PRINT  'Proaram  written  by:" 

600  PRINT  "LCDR  D.R,  Ferrell,  USN" 

610  PRINT  "Department  of  Aeronautics* 

620  PRINT  'Naval  Postgradiate  knool" 

630  PRINT  'Monterey,  California  93940  V«Septe«Der ,  1982)* 

640  DISP  ,,, 'Press  CONTINUE' 

650  PAUSE 

668  PAGESIZE  16  9  CLEAR 

670  !  Set  ice y  labels  for  composite  material  selection  display, 

630  ON  KEY*  l,'T300/5208*  GOTO  Ci  ■ 

690  ON  KEY*  2,'B(4)/5505'  GOTO  C2 

700  ON  KEY*  3. 'AS/3501*  GOTO  C3 

710  ON  KEY*  4,'SPiOQ2'  GOTO  C4 

720  ON  KEY*  5, "XV49/EP0XV  GOTO  C5 

730  ON  KEY*  6,"lher  Def'd"  GOTO  C6 

740  ON  KEY*  7.. 'References'  GOTO  441) 

750  ON  KEY*  3, 'SI  -)  ENG"  GOTO  349 

760  KEY  LABEL 

770  DISP  "  Selection  of  one  of  the  five  composites  listed  below  automatically* 

780  DISP  'reads  appropriate  engineering  constants  and  values  of  compressive,' 

790  DISP  'tensile,  and  snear  strengths  into  variable  names  for  later  comouta-* 

300  DISP  'tions,  All  values  are  in  31  units.  If  vou  anticipate  working  in" 

310  DISP  "English  units,  press  The  >SI  ->  ENG>  key.  If  vou  jesire  to  define' 

320  DISP  'a  material  not  listed  here,  select  'User  Def'd',  The  program  will" 

330  DISP  "then  interactively  reqsest  values  for  Ex,  Eu,  vx,  Es,  X,  <',  Y,' 

340  DISP  *Y',  and  3.  Select  'References'  for  an  explanation  of  the  orocnn' 

350  DISP  'and  a  listing  sf  references.' 

860  GOTO  360 

370  !  Set  data  pointer  to  read  appropriate  data  string  for  constants. 

380  Ci:  RESTORE  Dl  9  GOTO  RD 

390  C2:  RESTORE  D2  9  GOTO  RD 

909  C3:  RESTORE  D3  *  iJCTO  RD 

910  C4:  RESTORE  D4  9  GOTO  RD 
920  C5:  RESTORE  E  9  GOTO  RD 

°30  !  READ  md  DATA  statements  for  the  constants, 
949  RD:  READ  C3HPQSITE$,EX,EY,MUX.ES,X,XP/f,YP,S 
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950  IF  FLAG=0  THEN  GOTO  ii63 

960       EX=EX/SE  i  EY=EY/SE  9  E3=ES/SE 

970       X=X/5E  8  XP=XP/3E  9  Y=Y/SE  9  YP=YP/SE  9  S=S/SE 

980  GOTO  ii60 

990  !  DATA  statements,  See  oaaes  19  and  292,  ref  (a), 

1900  Di!  DATA  T300/S298, 131000000300, 10300003000,, 28,7173000000,1500800008, LSOIO 

03000,40000000,246008008,68000800 

1810  D2:  DATA  3(4)/SS0S, 204008033008, 18533303809, .23, 559)900009,1263099009,25300 

38303,61830088,202099390,67000008 

1320  D3;  DATA  AS/3501, 138030003830, 3969090000,, 3, 7100003800, 1447000080 ,144700030 

0,51703000,236000000,93000(380 

1030  D4:  DATA  SP1082, 33608000000,3279080000, .26,4143030000,1062303900,610000000, 

31003008,118000300,72080080 

1340  D5:  DATA  KV49/EP0XY, 76300030803, 5583008080, ,34,2300890000,1400008008,235000 

008,12008000,53003809,34003080 

1350  C6:  CLEAR 

i860  PRINT  "Input  name  of  the  Material:1 

1970  INPUT  COMPOSITE* 

1080  PRINT  9  PRINT  'Input  engineering  constant?: a 

1093  PRINT  "Ex,  Ey,  vx,  Es' 

1100  INPUT  EX,EY,NUX,ES 

iiiO  PRINT  USING  »2(3D . 3DE) , , 4D ,3» . 3DE"  ;  EX,EY,NUX,E3 

1120  PRINT  9  PRINT  "Input  compressive,  tensile,  and  shear  strengths:" 

1138  PRINT  "X,  V,   Y,  Y'  S" 

1140  INPUT  X,XP,Y,YP,S 

1158  PRINT  USING  '5<3D,3DE)"  ;  X,XP,Y,YP,3 

1160  CLEAR  9  PRINT  "Selected  composite  saterial  =>  ";COhPG3ITEl 

U73  PRINT  9  PRINT  9  PRINT  "Input  laminate  thidtness,  ho  (31  or  ENG)." 

1130  PRINT  8  PRINT  9  INPUT  hi 

1190  PRINT  8  PRINT  "ho  =";ho 

1208  WAIT  1800 

1210  CLEAR  9  DISP  9  DISP  'Do  you  wish  to  see  a  display  of  the  engineering" 

i229  DISP  "constants,  and  the  compressive,  tensile,  and  shear  contants?" 

i230  DISP  '   Y/N" 

1240  INPUT  A$ 

1253  IF  UPC*  (A*)="N"  THEN  GOTO  1420 

1263  CLEAR  >  pRINT  3  ?SINT  ■axni»««»«m«Bi«*«S«i««JW«tt««tt« 

mmmwimtmxunm*  *  print 

1270  PRINT  "Engineering  constants  and  compressive,  tensile,  and  shear11 

1238  PRINT  "constants  for  ";COHPOSITE$ 

1290  PRINT 

1300  PRINT  USING  "JX^A^IUDE"  j  ' 

1318  PRINT  USING  'i0X,SA,X,2D,3DE' 

1329  PRINT  USING  "iflX,SA,3X,,2D"  ;  ' 

1330  PRINT  USING  ■iOX,5A>2X>D.3DE" 
1340  PRINT  ?  PRINT  USING  "iiX,4A,4D, 
1350  PRINT  USING  "iOX.SA^D.SDE3  ; 
1360  PRINT  USING  9iiX .4A,2X, 2D . 3DE' 
1370  PRINT  USING  ■i0X)SA>X,3D.3DE'> 
1383  PRINT  USING  "11X,4A,2X,2D,3DE8 


Ex  =  ",EX 

JEy  =  \ET 

vx  =  \NUX 

"Es  =  ',ES 

3DE"  ;  "X  =  ' 

>x 

•X'  =  %XP 

.    <y  -   i  v 
>      '          ;• ' 

i   T  =  ',Yf 

;  *3  =  \S 
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1398  PRINT  §  PRINT  'u«HH«tn«u:n»au«<iH«tHw«mi«««iK««ti 
**»*«*«»**«*«*'  *  PRINT 

1409  DI5P  'Use  the  (ROLL)  xev  to  new  CRT  output.  ■,, "Press  CONTINUE* 

1410  PAUSE 

1420  CLEAR  9  DISP  i  MSP  ?  DI5P  "Coapotation  of  on-axis  modulus  and  coHpliance' 

1430  DISP  "terns.     Standby." 

1440  fi=i/(i-NUXA2*(EY/EX)) 

1450  QXX=M*EX  8  QYY=*ttO  ?  QXY=rt*MUX*EY  8  93S=EI 

1460  SXX=i/EX  ?  SYY=i/EY  8  SXY=-iNUX/EX)  3  3SS=i/ES 

1470  !  Compute  UiQ's,  the  linear  cowunations  of  on-axis  noduli  Oij,  ij=xvs! 

1480       UiQ=(3«QXX+3tgYY+2*QXYMIQSS)/8 

1490       U2Q=(QXX-QYY)/2 

1500       U30=<QXX+QYY-2*0XY-4«QSS)/8 

iSiO       U4Q=(QXX+QYY+««XY-4«SS)/8 

1S20       USQ=(QXX+9YY-2«QXY+4»Q5S)/8 

1530  !  Conpute  on-axis  A  Matrix  for  single  plj 

1540  AXX=QXX*hi  *  AYY=QYYtho  8  AXY=QXYtho  §  ASS=QSS*ho 

1550  !  Co«pute  inversion  of  A  natnx  (aij).     3ee  page  97, 

1560    5ET=ASS*(AXXtAYY-AXY12) 

1570  axx=AYY*ASS/0ET  9  a?y=AXXiASS/DET  9  axy=-<nXY*ASS/DET)  ?  iss=AXX*AYY/DET 

1580  DISP  9  DISP  'Do  <>n  wish  to  see  a  display  of  the  on-axis  8ij%  Sij's" 

1S9J  DISP  "(Jill's,  Aij's,  and  aij's?' 

1600  DISP  '   Y/N" 

1610  INPUT  A$ 

1&20  IF  UPC$  (A$)='N*  THEN  GOTO  1830 

1630  PAGESI2E  24  t  CLEAR 

1640  PRINT  i   PRINT  "ttt«t«»»««S»«««»«l«»»«*t«t«tt«t«<«$««*t»$ 

tmnmrntznnn'  ?  print 

1650  PRINT  JQn-axis  3ij;'s,  Sij's,  UiQ's,  Aij's,  and  nj's  for  '•.CCHPOSITE* 

1660  PRINT  '(Note;  Lj=xys)."  I  PRINT 

1670  PRINT  i  PRINT  USING  819X)2A>3D,2DE)2X>D,4DEl6X,A'  ;  '(  ',aXX,gXY, "H" 

1680  PRINT  USING  "5X, i8ft,iiX>2D.3DE,iX .A'  ;  'on-axis  Qij  =  C^YY,*!' 

1690  PRINT  USING  ,i9X,A,27X)D,4DE'  ;  ■(■tQSS 

1790  PRINT 

1710  PRINT  USING  "i9X,2A>D,4DE>2X,2D.3DE,6X)A"  j  '<  ",SXX>SXY,"Q" 

1720  PRINT  USING  "5Xyi8A,iiX>2D.3DE>6X>Aa  ;  'on-axis  Sij  =  t",3YY,T 

1730  PRINT  USING  "19X,A,27X>3D.2DE'  ;  '(',SSS 

1740  PRINT  8  SRINT  USING  "iSX^A ,2D ,  4DE"  ;  "U1Q  =  \M| 

1750  PRINT  USING  'iSX^A^D^DE'  ;  "U20  «  ',U2Q 

1760  PRINT  USING  'i3X,6A,21).4DE"  ;  *U3Q  =  ',U3Q 

1770  PRINT  USING  'i3X,6A,2D,4DE"  ;  'U4Q  =  \U4Q 

1780  PRINT  USING  ,13X)6A)25.4DE"  ;  'U5Q  =  *,U5Q 

1793  PRINT  ?  PRINT  USING  '19X>2A,2D.3DEf2X,.5DE16X)A*  ;  '(  VAXX.AXY.'Q' 

1300  PRINT  USING  '5X,i8A,iiX,D,4DE,6X,A"  ;  'on-axis  Aij  =  (",AYY,T 

1810  PRINT  USING  ,i?X)A,27X),5DE'  ;  VVASS 

1820  PRINT  9   PRINT  USING  "19X)2A,2D,3DE)2X,3D,2DE,6X,A*  ;  "<  " ;a<x,r*v , '?" 

1330  PRINT  USING  "SX , i 8A , i IX , 3D , 2DE , 6X , A*  ;  'on-axis  aij  =  (\ayy,"8" 

1840  3RINT  USING  -i9X , A ,27X ; 4D . DE"  ;  '(\ass 

1850  print  9  print  'mxitsmnmmmumimtmtmmmmtmumim 

tmmtntmnmi*  »  prist 

i860  DISP  'Use  (ROLL)  key  to  view  CRT  output.*,, 'Press  CONTINUE" 
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1871  PAUSE 

1388  CLEAR  9  PAGESIZE  16  9  DISP  9  DISP  9  DISP 

1390  DISP  'Computation  of  strenqth  parameters  in  strain  and  stress  space,' 

1900  DISP  'Stanabv.' 

i?10  FXX=i/(X*XP)  9  FYY=i/(YtYP)  9  FXY=FXYS*SQR  (FXXlFYY)  9  FSS=i/SA2 

1920  FX=i/X-l/XP  9  FY=i/Y-i/YP 

1930  GXX=FXX*QXXA2+2tFXY*8XX*0XY+fYY*9XYfc2 

1940  GYY=FXX*QXYA2+2*FXY*QXY*QYY+FYYt9YYA2 

1950  GXY=FXX*GXX*QXY+FXY*(QXX*QYY+9XY*2)+FYY$GXY*9YY 

1960  GSS=FSS*QSS*2 

1970  GX=FX*QXX+FY*QXY 

1930  GY=FX*QXY+FY*QYY 

1990  DISP  9  DISP 

2000  DISP  'Do  you  wish  to  see  a  display  of  strength  parameters  in  stress' 

2010  DISP  "space,  Fij,  and  strain  space,  Bij?" 

2020  DISP  '   T/N" 

2030  INPUT  A* 

2040  IF  UPC*  (A*)="N"  THEN  GOTO  2250 

2050  CLEAR  9  PAGESIZE  24  9  CLEAR 

2060  print  9  print  'mtmmmmtmnmttmmmtmmmttmmmm 
ttmnummmtv  i  print 

2070  PRINT  "Strength  parameters  in  stress  space,  Fij,  and  strain  space,  Bij." 

2D80  PRINT  'for  ";COHPDSITE$;*  (Note:     ij»xfs).'  9  PRINT  9  PRINT 

2090  PRINT  USING  "i7X12A,,5DE,2X,23.3DE,6X,Aa  ;  a(  a jFXXjFXY/O" 

2100  PRINT  USING  "tiX,7A,i4X>3D.2DE>6XJA"  ;  "Fij  =  C,FYY,"0" 

2110  PRINT  USING  'i7X,A>27X)3D,2DE8  ;  '(\FSS 

2120  PRINT 

2130  PRINT  USING  "12X,5A,2D,3DE'  ;  *Fx  =  ',FX 

2140  PRINT  USING  »i2X ,5A?2D.3I'E*  ;  "Fy  =  ",FY 

2150  PRINT  9  PRINT  i   PRINT 

2160  PRINT  USING  'i7X,2A,2D.3DE,2X,2D.3DE,6X,A'  ;  '<  '^XX^XY.'O' 

2170  PRINT  USING  "iiX^A, i4X,2D, 3DE,6X,A»  ;  "Gij  =  (*,GYY,'8' 

2180  PRINT  USING  'i7X,A,27X,2D,3DE'  ;  '(\GSS 

2190  PRINT 

2200  PRINT  USING  T2X,SA,D,4DEn  ;  *Gx  =  *,GX 

2210  PRINT  USING  "12X,5A,D.4DE"  ;  "Gy  =  ",GY 

2220  print  9  print  'mmmstmiimitmstmatmtimmmsmmttMm 
mmtnmmtmv  9  print 

2230  DISP  "Use  the  {ROLL}  key  to  new  CRT  output.', /Press  CONTINUE' 
2240  PAUSE 

2250  PAGESIZE  16  9  CLEAR  9  PRINT 

2260  3rint  9  print  'tmxiniinunxnmumutnmtmmmnnmmmn 

mmmtmnmna 

2270  print  *tttmnimmtnnitnt%mnmmmnnmnumnmmn%tn 

xmmmn*  9  print 

2280  PRINT  'The  program  will  now  compote  in-plane  and  flexural  properties  of" 
2290  PRINT  "laminates  specified  by  the  user,  Either  CORE  or  NO  CORE  Itminates' 
2300  PRINT  "can  oe  specified  and  any  number  of  plies  defined,  The  user  will* 
2310  PRINT  'note  that  this  program  takes  a  'brute  force'  approach  by  request-" 
2320  PRINT  'ing  an  entry  for  each  ply,  This  method  allows  simpler  program  " 
2330  PRINT  'code,  ensures  fewer  errors,  and  enhances  accuracy,' 
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2340  print  9  print  ntitmimmMtmtmmmtmtinnmttMtmmim 
mtmmnmmw 

2359  print  'tmximtminmmmtmmnmnxmmimittmmmmm 
mmmxtv  9  print 

2360  PRINT  'In-plane  coaptations  follow, ,." 

2370  PRINT  9  PRINT  "Will  the  lamnate  contain  CORE  or  NO  CORE?' 

2380  PRINT  "      C/HC  =>  END  LINE*  9  PRINT  9  PRINT 

2390  INPUT  A$ 

2403  IF  UPC$  (A$)='C'  THEN  GOTO  2SS0 

2410  CLEAR  *  PRINT  9  PRINT  JNG  CORE  lawinate:    Enter  >n>,  the  total  nunber' 

2420  PRINT  'of  plies,  surface  t»  surface,' 

2438  PRINT  u      n  =  ?" 

2440  INPUT  n 

2450  PRINT  »  PRINT  'n  =';n  9  PRINT  9  PRINT 

2460  h=n*ho  !  Thickness  if  NO  CORE  loiinatt. 

2470  PRINT  ?  PRINT  9  PRINT  "NO  CORE  lamnate  thickness,  h  ^";h 

2480  PRINT  9  PRINT  9  PRINT  "Njw  enter  the  ply  orientation    :  each  ply;" 

2490  PRINT  'THETA(i),  THETA(2),  ,,.   ,  THETA(n/2>" 

2500  PRINT  "XXX  --)  FROM  CENTER  TO  CUTER  SURFACE  <--  X&X    (IHPrjRT4HT ! ) " 

2510  PRINT  §  PRINT  §  PRINT 

2520  ti=t 

2530  GGSUB  THETA 

25^0  is=0  !  'Mune  fractisn  of  core, 

2550  GOTO  2300 

2560  CLEAR  *  PRINT  i  PRINT  "CORE  idMinate:    Enter  *n\  the  total  n«noer  of" 

2570  PRINT  "lamnate  plies,  and  'nc>,  the  total  nunoer  of  equivalent  core" 

2530  PRINT  "plies  that  define  its  thickness,  surface  to  surface,' 

2590  PRINT  "      n.nc  =  ?' 

2600  INPUT  n,nc 

2610  PRINT  9  PRINT  *n  =>;"      nc  =*;nc 

2629  rfAIT  1000  9  h=(n+ne)th«  !  ThicKness  of  CORE  laminate. 

2630  PRINT  i  PRINT  §  PRINT  'CORE  lawmate  thickness,  h  =J'.h 
2648  WIT  1000  9  CLEAR 

2650  PRINT  9  PRINT 

2660  PRINT  'Coontinq  outward  fro*  *he  center  of  the  laminate,  enter  an  ' 

2679  PRINT  "integer  nunber  ;'ti>  (t  initial)  that  reoresents  the  consecutive" 

2680  PRINT  'index  n outer  of  the  ply  adjacent  to  the  core,     (For  instance,' 
2690  PRINT  "if  the  core  is  a  total  of  4ho  thick,  count  2  equivalent  ply  layers" 
2700  PRINT  Jfro«  the  center  so  that  ti=3).' 

2719  PRINT  "      ti  =  ?' 

2729  INPUT  ti 

2739  PRINT  9  PRINT  9  PRINT  "ti  =";ti 

2749  mi  1090  9  PRINT 

2759  PRINT  *  PRINT  'Now  enter  the  pi?  orientation,  in  degrees,  for  each  plv:* 

2760  PRINT  'THETA(ti),  THETA(ti+l)   .,.   ,  THETA(n/2),' 

2770  PRINT  fttt  --)  CENTER  HOST  PLY  TO  OUTER  SURFACE  PLY  (-  *tt    (IMPORTANT) * 

2789  GOSUB  THETA 

2790  ;s=nctho/h  !    'olune  fraction  of  core, 

2890  CLEAR  *  MSP  9  CISP  9  DI3P  'Computation  of  in-plane  properties.    Standby," 
2810  !  CoMputation  of  in-plane  Modulus  of  lamnates. 
2329       ft(i,l)=tJiQX(n$ho)+Vi»U29+V2»U3Q 
2833       A(2,2)=Ui9»(n*h«)-ViWai+y2«J3Q 
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2840       A<i,2)=U4Qx(nXho>-V2xU3Q 

2350       A(6,6)=USQ*(n*ho)-V2XU3Q 

2860        A(i,6)=V3XU29/2+V4XU3Q 

2870       A(2,6)=V3XU2Q,'2-V4XU3Q 

2880  !  Conoutation  of  in-plane  compliance  of  laminates  (inversion  of  A ( i > j ) ) 

2890  DET=A(i)i)*A(2,2)«A(6>4)+2*A(i,2)tA^2,6;?A(i,6)-A(2>2)*A(i,6)A2-A(6)6):lcA(ll 

2)A2-A(i,i)tA(2)6)A2 

2900       a(i>i)=(A(2)2)tA(6)6)-A(2,6)A2)/DET 

2910       a(2,2)=(A(l,i)*A(6,6)-A(i,6)A2)/DET 

2920       a(i,2)=(A(i)6)<A(2>£>)-A(i>2)$A(S)t>))/DET 

2930       a(6J6)=(A(i>i)?A(2)2)-A(i,2)A2)/DET 

2940       a(i,6)=(A(i,2)tA(2,6)-A(2,2)tA(i,6))/DET 

2950       a(2,6)=(A(i,2)*A(i,6)-A(i,i)tA(2,6))/DET 

2960  I  Computation  of  normalized  in-plane  Modulus  and  compliance 

2970  FOR  1=1  TO  6 

2980       FOR  J=l  TO  6 

2990  AS(I,J)=A(I,J)/i 

3000  as(I,J)=a<I,J)xh 

3010       NEXT  J 

3020  NEXT  I 

3830  !  Computation  of  effective  in-plane  engineering  constants. 

3040       Eio=i/(a(i,i)Xh)  !  Longitudinal  modulus 

3050       E2o=i/U(2,2)*h)  !  Transverse  modulus 

3060       E6o=i/\a(6,6)Xh)  !  Shear  modulos 

3070       v2i«=-(a<i>2)/a<i>i>)  (  Poisson's  ratio 

3080  DISP  ?  DI3P  i  DISP  "Do  sou  wish  to  see  a  display  of  effective  in-plane  " 

3090  DISP  'engineering  constants.,  the  in-plane  stiffness  matrix  (Aij)" 

3103  DISP  'and  its  inverse  <aij),  and  the  normalized  versuns  of  these' 

3110  DISP  'fwtrices  (Ai j*  and  aij*)?" 

3120  DISP  *     Y/H' 

3130  INPUT  A* 

3140  IF  UPC*  (A$)='N'  THEN  GOTO  3430 

3150  PAGESI2E  24  9  CLEAR  ?  PRINT 

3160  print  §  print  attx<txxxm*xmx*xxxxmxxxxmm*m»x?*x*mm*3xxixtx 
tmmmtmttmv  ?  primt 

3170  PRINT  "In-olane  parameters:    Aij's,  aij's,  enaineering  constants,  and" 
3130  PRINT  "normalized  Aijx's  and  aijX's  for  ■; COMPOSITE* ;"     (Note:     ij=i26)' 
3190  PRINT  i  PRINT 

3200  PRINT  USING  "i9X.,2A,,2D,3DE,2X,D,4DE,2X)D,4DE"  ;  '(  * , A< i , i) , A< i >2) , A(i ,6) 
3210  PRINT  USING  a<X,i7A,i3X,2D,3DE,2X,D,4DE"  J  "in-plane  Aij  =  r,A<2,2),A(2,e! 

3220  PRINT  USING  'M^^D^DE'  ;  'C,A(6,i) 

3230  PRINT 

3240  PRINT  USING  "i9X,2A>2D.3DE,2X,D.4DE,2X,B.4DE"  \  '(  ,,a(i)i)>a(i,2)>a(i,6) 

3250  PRINT  USING  '4X,i7Afi3X)2D,35E,2X>D.4!)E'  ;  "in-plane  aij  =  (\a(2,2),a(2,6i 

3260  PRINT  USING  *i?Xf A >27X,D . 4DE"  ;  ,C>a(6,&) 

3270  PRINT 

3280  PRINT  USING  ,13X,oA,D.4DE'  ;  "Eio  =  \Eii 
3290  PRINT  USING  "13X^,5, 4DE"  ;  "E2o  =  \E2o 
3300  PRINT  USING  "i2X,7A,.Sl>"  ;  "v2io  =  ",»21o 
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3310  PRINT  USING  'ilX.&A.D.JDE"  ;  "E6o  =  ",E6o 

3320  PRINT 

3330  PRINT  USING  'i?X,2A,D.4DE,2X,i>.4DE,2X,D.4DEtt  ;  '<  ,,AS(i)i)JAS(i,2),AS(i,&) 

3340  PRINT  USING  "3X,i7A,i4X>D.4DE,2X>». 4DE"  ;  a  in-plane  Aij*  =  C , AS( 2,2) ,AS( 2 , 
S) 

3350  PRINT  USING  "i9X>A}27X ,D. 4DE"  ;  '(\AS(6,6) 

3360  PRINT 

3370  PRINT  USING  ,i?X)2A)D.4DEi2X>D.4DE)2X,D,4DE"  ;  u(  ,,05(i,i))as(i)2),os(i>6) 

3330  PRINT  USING  a3X,i7A,i4X,D.4DE,2X1D.4DE'  ;  'in-plane  aij*  =  (,,«(2,2),as<2l 

6) 

3390  PRINT  USING  ,i9X>A)27X)D.4DEB  ;  ,(,,os(6  6) 

3400  print  %  print  *xxxxxxxxxxx%ix%t%xxxxxxxxxxxxxxxx%x%xxxx%xx%itx%tx%x%%xtxtx* 
xxtxxxxxxxxnxxxxxxx*  i  print 

3410  DISP  'Use  (ROLL)  key  to  view  CRT  output. ',,'Press  CONTINUE' 

3420  PAUSE 

3430  CLEAR  8  PAGESIZE  16 

3448  DISP  'It  will  now  De  necessary  to  input  the  loads,  Ni,  N2,  and  N6." 

3450  DISP  'Consider  two  options:* 

3460  DISP  '       Case  (a);     Input  selective  unit  loads  to  determine' 

3470  DISP  *  naxiflun  loading  allowaoie." 

34S0  DISP  "       Case  (b):     Input  actual  loading  to  determine  strain" 

349!)  DISP  '  invariants  and  strength  ratios.1 

3S00  DISP  '(The  progran  operates  the  sane  m  either  case.     It  is  up  to  the  ' 

3510  DISP  "user  to  interpret  the  output  properly).' 

3523  PRINT  §  PRINT  9  PRINT  'Input  loads  Ni,N2,N6"  9  PRINT  9  PRINT 

3S30  INPUT  Ni,N2,N6 

3S40  PRINT  9  PRINT  'Ni  =';Ni;*      M2  =';N2;'      N6  =B;N6 

3559  WAIT  1301 

3560  !  Compute  the  in-plane  strain  variables,  eio,  e2o,  e6o,  (T4.2,  p . 121) 
3S70       eio=a(i;i)INi+aa>2)*N2+a(i>6/*H6 

3S30       e2o=a(i,2)*Ni+a(2,2)«!2+<i<2,6)»N6 

3599  efco=a(i,6)*Ni+a(2,6)*N2+a(6,S)*N6 

3600  PRINT  »  PRINT  9  PRINT  I  CLEAR 

3610  PRINT  %  PRINT  'Now,  enter  THETAt,  ^he  orientation  if  the  pl»" 
3620  PRINT  "to  be  exanined  (in  degrees).' 
3630  PRINT  '   THETAt  =  ?■  \   ?RIMT  *  PRINT 

3649  INPUT  THETAd 

3650  PRINT  9  mm  'THETAt  =';THETAd 
3668  JAIT  1000 

3670  THETAt=THETAd*PI  /ISO 

3633  CLEAR  9  DISP  9  DISP  §  DISP  'Computation  of  in-plane  properties.    Standby" 

3690  !  Compute  the  on-axis  strain  variable,  exo,  eyo,  eso,  for  the  given  ply. 

3700  '  T2.6,  p. 53 

3713       p=(eio+e2o)/2  9  g=(eio-e2o)/2  9  r=e6o/2 

3720       *«=COS  (2ITHETAt)  ?  nn=SIN  (2tTHETAt) 

3730  exosp+qtiw+rtnn 

3740  e«o=p-gt««-rtnn 

3753  es»=-(2tq<nn)+2tr*«n 

3760  !  Compute  the  on-axis  stress  variables,  sxo,  syo,  sso,  for  the  given  ply, 
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3770  '  Ti.6,  p. 14 

3780  5xo=QXX*exo+QXYtevo 

3790  sy»=QXYiexo+9YY*eyo 

3800  55O=QSS*€50 

3810  GQSUB  SSSR     !  Go  subr  to  conpute  the  stress-strain  strength  ratios,  R&R» 

3820  GQSUB  5UAD    !  Go  subr  to  conpute  the  quadratic  strength  ratios. 

3830       sto=Rq/h  S  stopp=Rqp/h 

3840  CLEAR  i  PAGESIZE  24  i  CLEAR  ?  PRINT 

3850  print  i  print  'txtxtxxxtxtxttxtxixxxxxxxxxxtxtxxxxxxxxxxxxtxtxxxxxxxxxxxxxx 
txxxxxxxxxxxxmmx*  §  print 

3360  PRINT  USING  U10X,23A)2D.3DE1  ;  'in-plane  strain:  elo  =  ",eia 

3370  PRINT  USING  "27X...SA.D .4DE"  ;  *e2o  =  *,e2o 

3880  PRINT  USING  "27X,6A,D.4DE"  ;  "e6o  =  ',e6o 

3890  PRINT  *  PRINT  'For  THETAt  =';THETAd;'  ....  (»; COMPOSITE* ;')"  ?  PRINT 

3900  PRINT  USING  ,20X>22A)D.4DE'  ;  'on-axis  strain:  exo  =  ",exo 

3910  PRINT  USING  a36Xf6A>29.3DE"  ;  "eyo  =  \eyo 

3920  PRINT  USING  *36X,6A,D.4DE"  ;  *eso  =  ",eso 

3938  PRINT 

3940  PRINT  USING  '28X,22A,2D.3DEB  ;  "on-axis  stress:  sxo  =  \sxo 

39S0  PRINT  USING  'SiX^A^D. 3DE"  ;  "syo  =  ',syo 

3960  PRINT  USING  >36X>6A,D.4DE*  }  "sso  =  \sso 

3970  PRINT  3  PRINT 

3980  PRINT  USING  'lOX^SA^X^iiX^A'  ;  'strength  ratios","*'/*" 

3990  PRINT 

4000  PRINT  USING  ,15X,10A,4X,D,DDE)3X,D,DDE'  ;  "*ax  strain  ',Re,Rep 

4010  PRINT  USING  ,1SX,10A,4X)D,:)DE)3X,D.DDE"  ;  **ax  stress  ",Rs,Rsp 

4020  PRINT  USING  "iSX> iOA,4X>D . DDE.3X.D . DOE*  ;  "quadratic    ",Rq,Rqp 

4030  °RINT 

4040  PRINT  USING  '20X,24A,1).4DE'  ;  "averaqe  stresses:  sto  =  Vsto 

4050  PRINT  USING  '37X^,0. 4DE'  ;  'sto'  =  \stopp 

4060  PRINT  USING  "30X,22A"  ;  '<=f(quadratic  R  i  R>))" 

4079  print  9  print  'xxxxxxtxtxxtxxxxxtxxxtxxxxxxxixmxxxxxxxxxxixxmxtxxxnm 
mtttmmmtmt'  t  print 

4080  MSP  'Use  (ROLL)  key  to  new  CRT  output.' ,, 'Press  CONTINUE' 
4090  PAUSE 

4180  CLEAR  §  PAGESIZE  15 

4110  DISP  'Do  you  wish  to  conpute  in-plane  properties  for" 

4120  DISP  'a  new  load  and  THETAt?' 

4130  DISP  '   Y/H" 

4140  INPUT  A* 

4150  IF  UPC*  (A*)='N'  THEN  GOTO  4201 

4160  DISP  *  DIS?  "Enter  new  Hi/fC,^.' 

4170  INPUT  Ni,N2,N6 

4180  PRINT  ?  PRINT  2  PRINT  'New  loaa:  Nl  =',N1;'   N2  =";N2;'   N6  =';N6 

4190  WAIT  iOOO  ?  GOTO  3568 

4200  CLEAR  §  PRINT 

42io  print  §  print  •%xtxxxtxmxttx%xxxxxx%xxxtxixxxixxxxxxxxxxxxtx%xxxxximxxx 

tmmttmmxtttf 

4220  print  'ttmmmmttttttmtmttmxttmmmttmmimttmmmt 

ttmtmuv  ?  print 

4230  PRINT  'The  program  will  now  cowoitte  flexural  priperties  for  the  laninate" 
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4240  PRINT  "defined  during  the  in-plane  definitions.    Properties  calculated" 
4250  PRINT  "include  the  flexural  stiffness  «atnx  (Dip,  the  flexural  comoli-" 
4260  PRINT  "ance  matrix  (dij),  effective  flexural  hamate  moduli  (Eif),  and" 
4270  PRINT  "the  normalized  versions  of  stiffness  matrices  (Dij*  and  dij*)." 

4290  print  e  print  'mtmummmmsttmmimtmsttmtmtmmmtt 
mmtmmtmtw 

4290  PR  INT  *  mmttmmmtmimmmttmmmnmtmummmmm 
************ * 

4300  DISP  i   DISP  ,,, "Press  CONTINUE" 

4310  PAUSE 

4320  CLEAR  §  DISP  3  DISP  %  DISP 

4330  DISP  'Comoutation  of  flexural  properties.  Standby." 

4340  hs='i-zsA3)*hA3/i2 

4350  !  Compute  flexural  Modulus  (T5.3,  p. 181) 

4360       D(i>i)=lliQ*hs+ViD*U2Q+Y2Dt038 

4370       D(2,2)=UiQ*hs-VlDIU2Q+V2D»U38 

4380       D<i,2)=U4Q*hs-V2D*U3Q 

4390       D(6,6)=US9*hs-V2D*U3Q 

4400       D(i,6)=U29*V3D/2+V4D*U3Q 

4418       D(2,6)=U2Q*V3D/2-y4D*U3Q 

4*29  !  Compute  flexural  compliance  (inverse  of  Dij) 

4430       DET=D(i,i)*D(2)25*D(6)(»)+2*D(i)2)*D(2)6)iD(i)6)-D(2)2)*D(l>i>)A2-D(6,6)*D 

(i,2)*2-D<i,i)t5(2,6)A2 

4440   d(i,l)=(D(2>2)*D(6)6)-D(2)6)A2)/DET 
4450   d(2,2)=(D(i,i)*D(6,6)-D(i,6)A2)/DET 
4460   d(i,2)=(D(i>6)*D(2>6)-D(i,2)*D(6)6))/DET 
4470   d(6)6)=(D(i>i)*D(2)2)-D<l,2)A2)/DET 
4430   d(i,6)=(D(i,2)*D(2,6)-D(2>2)*D(l)6))/DET 
4490   d(2,6)=(D(i,2)tt(i,6)-D(i,i)*D(2,6))/DET 
4500  !  Compute  normalized  versions  of  Dij  and  dij. 
4510  FOR  1=1  TO  6 

4529  FOR  J=i  TO  I 

4530  DS(I,J)sD(I,J)*i2/h*3 
4548  ds(I,J)=d(I,J)*hA3/i2 
4559   NEXT  J 

4S68  NEXT  I 

4570  !  Compute  effective  flexural  laminate  moduli. 

4530   Eif=i2/(hA3*d(i>l)) 

4593   E2f=i2/(hA3*d<2,2)> 

4600   E6f=i2y(hA3*d(6,6):' 

^610   »2if=-(d<i)2)/d(i,i)) 

4623  DISP  I   DISP 

4630  DI'SP  'Do  viu  wish  to  see  a  display  of  flexural  modulus  and  compliance* 

4640  DISP  "matrices  iDij  and  dij),  effective  flexural  laminate  moduli,,  lnd" 

4650  DISP  'th*  normalized  versions  of  the  modulus  and  compliance  matrices)" 

4660  DISP  '   Y/N" 

4670  INPUT  A$ 

i&SO  IF  UPC$  iA*)='N"  THEN  GOTO  4970 

4690  PAGESIZE  24  t  CLEAR  ?  PRINT 

4708  PRINT  i  PRINT  't******************************************************!?*** 

**?***?**«*********"  t   PRINT 
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4713  PRINT  "Flexural  properties  Dij,  dij,  affective  flexural  engineering" 

4720  PRINT  "constants,  and  the  normalized  Dij*  and  dij*.  (Note:  ij=i26)' 

4730  PRINT 

4740  PRINT  USING  "19X > 2A ,D . 4DE V2X, . SDE,2X ,D ,4DE"  ;  *<  ",0(1,1), D(i,2),D(l,6) 

4750  PRINT  USING  "4X,16A,14X,D.4DE,2X,D.4DE"  ;  "flexural  Dij  =  (",D(2,2),D<2,6) 

4760  PRINT  USING  '19X,A,27X,D.4DE'  ;  n(',D(6,6> 

4770  PRINT 

4780  PRINT  USING  'i9X,2A,D.4DE,2X)D.4DEi2X,D.4DE*  j  "(  ,,d(i)i),d(i)2))d(i,6) 

4798  PRINT  USING  "4X,i6A,i4X,D,4DE,2X,D.4DE"  ;  "flexural  dij  =  <\d(2,2),d(2,6i 

4800  PRINT  USING  '19X,A,27X,D,4DE"  ;  '(",d(6,6) 

4318  PRINT 

4829  PRINT  USING  "i3X,6A>D.4DE''  j  "Eif  =  ",Eif 
4833  PRINT  USING  "13X,7A,D,4DE"  ;  'E2f  =  ",E2f 
4840  PRINT  USING  "i2X,9A,.5D"  ;  "v2if  =  ",v2if 
A8SQ  PRINT  USING  '13X,6A,D.4SEB  ;  "E6f  =  ",E6f 
4860  PRINT 

4379  PRINT  USING  "i9X,2A,D.4DE,2X,D.4DE,2X,D.4DE"  j  "(  ",DS(i,i),DS(i,2),DS(i,6) 

4830  PRINT  USING  "3X,17A,14X,D.4DE,2X,D.4DE"  ;  "flexural  Dij*  =  (VDS(2,2),DS(2, 
6) 

4390  PRINT  USING  "19X,A,27X,D.4DE"  ;  '(",D3(6,6) 

4900  PRINT 

4913  PRINT  USING  'i9X,2A,D,4DE,2X,D,4DE,2X,D.4DE'  ;  '(  ,,ds(i,i).,ds(i>2),ds(l,6) 

4920  PRINT  USING  "3X,17A,i4X,i).4DE,2X,D.4DE"  ;  "flexural  dij*  =  (,,ds(2,2),ds(2, 

6) 

4933  PRINT  USING  ,i9X)A,27X,D,4DE'  ;  "<",ds(6,6) 

4940  print  §  print  'tmmitmmmmtxnttimmmmwmtnxmmtm 
nmmmmumr  i  print 

4958  DISP  "Use  (ROLL)  key  to  new  CRT  output,",, "Press  CONTINUE' 

4960  PAUSE 

*979  PAGE3IZE  16  8  CLEAR 

4?89  DISP  'It  will  now  be  necessary  to  input  nonents  Hi,  H2,  and  H6," 

4990  DISP  "Consider  two  options.' 

5000  DISP  '   Case  (a):  Input  selective  unit  wonents  to  determine  the* 

5310  DISP  '  «ax  loading  allowable." 

5920  DISP  '   Case  (b):  Input  actual  nowents  to  determine  strength  ratios' 

5033  DISP  "  and  curvatures." 

5940  DISP  "(The  program  operates  the  sane  in  either  case.  It  is  uo  to  the' 

SQ50  DISP  "user  to  property  interoret  the  results).* 

5060  PRINT  9  PRINT  'Input  undents  Ht,M2,H6"  %  PRINT  «  PRINT 

5073  INPUT  H1,M2,H6 

3030  PRINT  ?  PRINT  'Mi  =";Hij"   H2  =";h'2;*   H6  =";rt6 

5099  WAIT  1933 

5i00  CLEAR  i  PRINT  8  PRINT 

5ii9  PRINT  'Now,  enter  the  index  number,  t,  and  the  orientation,  THETAt," 

5120  PRINT  'of  the  ply  to  be  examined,' 

5130  PRINT  '   t,  THETAt  =  J"  %  PRINT  g  PRINT  §  PRINT 

5143  INPUT  tJHETAd 

5150  PRINT  »  PRINT  't  =■;*;"   THETAt  =";THETAd 

5163  WAIT  1303 
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5170  THETAt=THETAd*PI  /180 

5130  CLEAR  9  DISP  g  51SP  8  DISP  '  Standby1 

5190  !  Computation  of  curvatures,  ki,     (T5. 2,  p. 1761 

5200       k l=d  <  i , i ) *Mi+d ( i  s  2 ) *M2+d  <  i  ?  S ) *M6 

5210       k2=d(i,2)IHi+d(2,2)tM2+d(2,6)tt1> 

5223       k6=d(i,i)tti+d(2,6)*M2+d(6;i)«H6 

5230  eio=kiit*no  !  actually  elf;  called  eio  for  sobr  QUA] 

5240  e2o=k2*t*ho 

5250  e6o=k6tt*ho 

5268  eif=€io 

5270  e2f=e2o 

5280  e6f=e6o 

5298  !  Compute  jff-axis  compliance  for  stress  calculations. 

5308       Vi=hotCQS  (2$THETAt) 

5319       V2=h«*COS  <4*THETAt) 

5329  V3=ho«SM  (2»THETAt) 

5330  U4=hotSIM  (4«THETAt! 

5340  A(l,i)=UiQ*ho+Vi*U2Q+V2*U3Q 

5359  Q<l,l)=A(i,i)/ho 

5369  A(2,2)=UiQ*ho-Vi*U29+V2*U3Q 

5370  Q(2,2)=A(2,2)/h« 
5380            A(iv2)=lJ4Qiho-V2tU3l) 
5390  Q(i,2)=AU,2)/h» 
5400            A(6,6)=U5Q*ho-V2*U3< 
5410                 Q(6)6)=A(6,6)/ho 
5420            A(i)6)=U2QtV3/2+y4tU3l 
5430                 Q(i,6)=A(i>6)/ho 

5449  A(2)6)=028»V3/2-V4*U3li 

5450  9(2)6)=A(2)6)/ho 

5469  sif=Q(l)l)Jeif+Q(i>2)»e2f+Qa)6)Ie6i 

5470  52f=9(i,2)teif+Q<2,2)*e2f+Q(2,i)*e6f 
5480  s6f=Q(l,6)te6f+Q(2,6)le2f+9(6j6)*e6f 
5490  !  Cortpute  on-axis  strain  for  R/R'  calculations. 
5500  !   (T2.6,  p. 53) 

5510       p=(eio+e2»)/2  9  q=(eio-e2o)/2  8  r=et>o/2 

5529  mm=COS  (2tTHE7At)  §  nn=S!N  (2*THETAt) 

5530  exo=p+q*n«+r*nn 
5540  eyo=p-q*mn-r*nn 
5550            eso=-(2tq*nn)+2*r*wi 

5560  !  Compute  on-axis  stress  for  R/R*  calculations. 
5570  sx»=QXX*exo+9XY*eyo 

5589  3?»=QXYtex«+flYY*eM 

5590  5so=QSSteso 
5600  GOSUB  SSSR 

5610  GOSUB  QUAD 

5620        sto=6ttq/V2  §  stopp=6*Rqp/hA2 

5630  CLEAR  9  PAGESIZE  24  9  CLEAR 

5640  PRINT  i  PRINT  *ntttttttttmwtittunmnmt%ttm%nutunututitin 

mmmmmmxi*  i  print 

5659  PRINT  USING  "10X . 16A.D . 4DE*  ;  'curvature ,  ki  =  Mti 

5660  PRINT  USING  "2iX,5A,9.4D£"  :  "k2  =  *,*2 
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;  "flex 

'1  strains, 

eif  = 

■>«t» 

;  "e2f  - 

'■   *,e2f 

j  "e6f 

=  \2&f 

;  "flex 

4  stresses 

,  sif  = 

■>sif 

=  ',s2f 

;  "s6f  : 

■■   ',s6f 

A1  .  ■«■ 

hrann  th  */*  +  i 

5670  PRINT  USING  "2iX,5A,D.4DE"  ;  "k6  =  ",k6 

5683  PRINT  9  PRINT  "Fop  t  ="jt}*and  THETAt  ='jMTftd  I  PRINT 

S690  PRINT  USING  "4X ,22A,D. 4DE" 

5700  PRINT  USING  *23X)6A)D.4DEa 

5710  PRINT  USING  l2IX1MlMKB 

5720  PRINT 

5730  PRINT  USING  •3X,23A,D.40E" 

5740  PRINT  USING  ,28X,6A,D.4DE' 

S750  PRINT  USING  B20X,6A,D.4DE* 

5760  PRINT 

S770  PRINT  USING  " ISA f3X, A,i iX,2A*  ;  "strenqth  nt iosV -R-.'R - ' 

5738  PRIST 

5790  PRINT  USING  "SX,i3A,4X,D.DDE,3X,D,DDE"  ;  "wax  strain  ',Re,Rep 

5800  PRINT  USING  ■SX,idA>4X,D.DDE,3X,D.DDE"  ;  "«ax  stress  ',Rs,Rsp 

5810  PRINT  USING  "SXf iOA^XjD.DDE^X,!) .DDE"  ;  'quadratic    \Rq,Rqp 

5320  PRINT 

5830  PRINT  USING  "20X>24A)D.4DE"  ;  "average  stresses:  sto  =  J,sto 

5840  PRINT  USING  "37X,7A,D.4DE"  }  "sto'  =  ",stopp 

5850  PRINT  USING  "30X,22A"  ;  "(=f< quadratic  R  4  R'))" 

5860  print  9  print  •tmtnmmmtmnmmmtmmmtmmtmmtim 
tummttttmutf  §  print 

5879  DISP  "Use  (ROLL)  key  to  new  CRT  output, './Press  CONTINUE" 

5380  PAUSE 

5390  CLEAR  8  PAGESIZE  16  9  DISP  §  DISP  §  DISP 

5900  DISP  "Do  you  wish  to  compute  a  new  set  of  nonents,  Hi,H2,H6,  index  t," 

5910  DISP  "and  orientation,  THETAt?  (NO  will  terninate  proqran  operation)." 

5920  DISP  '   Y/H" 

5930  INPUT  A$ 

5?40  IF  UPC*  (A$)="N"  THEN  GOTO  6000 

5950  PRINT  ?  PRINT  9  PRINT  "Input  new  loads  Mi,H2,H6"  ?  PRINT  9  PRINT 

5960  INPUT  HI, (12 ,,16 

5970  PRINT  8  PRINT  "New  aorcnts:    Hi  =";Hi;'      H2  =";H2;"      H6  =';H6 

5980  WAIT  1090 

5990  GOTO  5100 

bOOO  PRINT  9  PRINT  "End  of  jrograw.     Thank  |iott.' 

6910  ENS 

632?  ! 

6039  ! 

6040  THETA:  !  Subroutine  to  enter  THETA's  for  calculation  of  Ws  and  ViD's. 
6953  FOR  t=ti  TO  n/2+ti-i 

6860       PRINT  9  PRINT  "      THETAC;  t- a)  =  ?«  9  PRINT 

6079  INPUT  THETA(t) 

6080  PRINT  ,THETA(t) 

6093       THETAt t)=THETA(t)*PI  /IB  J 

Si 00        GOSUB  Vi 

6110       GOSUB  WiP 

6120  NEXT  t 

6i30  x=2tht< 

ii49  Vi=x*Vi  9  V2=x*V2  9  ^2=xtV3  9  V*=xlV4 

6150  x=2/3lho*3 

al60  U1D=t*V1D  9  V2D=xtV2D  9  V3D=xtV3B  9  V4D=x*V4D 

6179  RETURN 
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6188  ! 

6190   Vi: 

6290     Vi=Vl+COS  (2*THETA(t)) 

6219  V2=V2+C0S  <4*THETA(t)) 

6220  V3=V3+SIN  (2ITHETA(t)) 
6230     V4=V4+SIN  (4*THETA(t)) 

6249  RETURN 

6250  ! 
6260   ViD: 

6279  x=tA3-(t-i)A} 

6280  ViD=VlD+COS  (2tTHETA(t))*x 
6290  V2D=V2D+C0S  (4lTHETA(t))*x 
6300  V3D=*/3D+SIN  (2*THETA(t))*x 
6310  V4D=V4D+SIN  (4lTHETAit) )?x 
6328       RETURN 

6330  ! 

6340  ' 

b350  QUAD:   !  Transforation  of  quadratic  strength  parameters  in  strain  space 

6360  !  in  Multiple  angle  functions.    T7.6a,  p .298 

6375       mh=C0S  (THETAt)  »  nn=SIN  (THETAt) 

6339  Gii=GXX*MH*4+GKY«nn*4+GXY*2tiwA2«nn*2+GS&l4ttt«A2tnn12 

6398  G2^GXXtnnA4+GYY*w^4+GXv*21^•^«A2tnnA2+GSSt4t«MA2^nnA2 

6400  Gi2=GXXt«MA2tnnA2+GYY»««A2JnnA2+GXYt(«nA4+nnA4)-G3Si4*nNA2l£nnA2 

6410  G66=5XX«H«A2tnnA2+GYY?M«A2*nnA2-GXYt2?««A2lnnA2+GSS*(rt«A2-nnA2)A2 

6420  G16=GXXJfwA3*nn-GYY*HfttnnA3+GXY*(««innA3-^A3tnn)+G33*2t(i«*nnA3-flftA3Jnn) 

5439  G26=GXX««««nnA3-GYY*«nA3*nn+GXY»(««A3*nn-fifl*nnA3)+GSS*2*(«MA3tnn-rt«lnnA3) 

6449  Si=GX*««A2+GY*nn*2  !  17.6b,  p .298 

•6459  G2=GX*nnA2+GY*m-A2 

6460  G6=GX*«nlnn-GTtMfllns 

6470  a=GilleioA2+G22?e2oA2+2»Gi2iela*e2o+G66te6oA2+2tG16<eio<ebo+2<G26<e2o*e6o 

6480  Q=Glteio+G2*e2oiG6*e6« 

6499  c=-i 

^500  Rq=(-ts+S88  <bA2-<>»qfc>}/(2ta) 

6510  Rap=ABS  ((-b-SSR  '.bA2-4*a*c))/(2*a)) 

6529  RETURN 

6530  ' 
6540  ! 

bSSfl  SSSR:   !  Subroutine  to  conoste  wax  stress-strain  R  ana  R' 

6569  Res=S/AB£  (esolEB) 

6579  IF  exo>()  THEN  Rex=X/(exo*EX)  ELSE  Rex=<P/ABS  (exo*EX) 

6380      IF  ?'Jo>9  THEN  Rey=Y/(eyo*EY)  ELSE  Rey=YP/ABS  <evo*EYi 

6590         IF  Rex (Rev  THEN  Re=Rex  ELSE  Re=Rey 

6609  IF  ResCRe  THEN  Re=Res 

6619  IF  ex»<9  THEN  Rex=X/ABS  (exo«EX)  ELSE  Rex=XP/(exolEX) 

6620  IF  eyo<9  THEN  Rey=Y/ABS  leyolEY)  ELSE  Rey=YP/(eyo*EYi 
6o30  IF  Rex<Rey  THEN  Rep=Rex  ELSE  Rep=Rey 

6648  IF  Rss<Rep  THEN  Reo=Res 

6650  Rss=S7ABS  (sso) 

6660  IF  5xo>0  THEN  Rsx=X/sxo  ELSE  Rsx=XP/ABS  (sxo) 

6670      IF  syo>0  THEN  Rsy=Y/syo  ELSE  Rsy=YP/ABS  (syo! 

6680         IF  Rsx^Rsy  THEN  Rs=Rsx  ELSE  Rs=Rsy 
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S690  IF  Rss<R5  THEN  Rs=Rss 

6700      IF  sxo(0  THEN  Rsx=X/ABS  (sxo)  ELSE  Rsx=XP/sxo 

6710         IF  syo<0  THEN  Rsy=Y/ABS  (svi)  ELSE  Rsv=YP/svi 

6720  IF  Rsx<Rsy  THEN  Rsp=Rsx  ELSE  Rsp=Rsy 

6730  IF  Rss<Rsp  THFJ  Rsp=Rss 

6740  RETURN 

67S0  END 
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